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I. INTRODUCTION 



Much of the current understanding of solid state physics is based on a picture of non- 
interacting electrons. This is clearly true at the elementary level, but in fact extends to many 
areas of current research, examples being the physics of disordered systems or mesoscopic 
physics. The most outstanding examples where the non-interacting electron picture fails 
are provided by electronic phase transitions like superconductivity or magnetism. However, 
more generally one clearly has to understand why the non-interacting approximation is 
successful, for example in understanding the physics of metals where one has a rather dense 
gas (or liquid) of electrons which certainly interact via their mutual Coulomb repulsion. A 
first answer is provided by Landau's theory of Fermi liquids |l]-[3| which, starting from the 
(reasonable but theoretically unproven) hypothesis of the existence of quasiparticles shows 
that the properties of an interacting system of fermions are qualitatively similar to that 
of a non-interacting system. A brief outline of Landau's theory in its most elementary 
aspects will be given in the following section, and a re-interpretation as a fixed point of a 
renormalization group will be discussed in sec.p|. 

A natural question to ask is whether Fermi liquid like behavior is universal in many- 
electron systems. The by far best studied example showing that this is not the case is the 
one-dimensional interacting electron gas. Starting with the early work of Mattis and Lieb 
(U, of Bychkov et al. |J and of others is has become quite that in one dimension Landau 
type quasiparticles do not exist. The unusual one-dimensional behavior has now received 
the name "Luttinger liquid". This is still a very active area of research, and the rest of 
these notes is devoted to the discussion of various aspects of the physics of one-dimensional 
interacting fermions. 

The initial plan for these lectures was considerably wider in scope. It was in particular 
considered to include a discussion of the Kondo effect and its non-Fermi-liquid derivatives, 
as well as possibly current theories of strongly correlated fermions in dimension larger than 
one. This plan however turned out to be overly ambitious and it was decided to limit 
the scope to the current subjects, allowing for sufficiently detailed lectures. Beyond this 
limitation on the scope of the lectures, length restrictions on the lecture notes imposed 
further cuts. In view of the fact that there is a considerable and easily accessible literature 
on Fermi liquid theory, it seemed best to remain at a rather elementary level at this point 
and to retain sufficient space for the discussion of the more unusual one-dimensional case. 
It is hoped that the references, especially in the next and the last section will allow the 
interested reader to find sources for further study. 



Landau's Fermi liquid theory is concerned with the properties of a many-fermion sys- 
tem at low temperatures (much lower than the Fermi energy) in the normal state, i.e. in 
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the absence of or at least at temperatures above any symmetry breaking phase transition 
(superconducting, magnetic, or otherwise). The ideal example for Landau's theory is liquid 
helium 3, above its superffuid phase transition, however, the conceptual basis of Landau's 
theory is equally applicable to a variety of other systems, in particular electrons in metals. 
Quantitative applications are however more difficult because of a variety of complications, 
in particular the absence of translational and rotational invariance and the presence of 
electron-phonon interactions, which are not directly taken into account in Landau's theory. 
Subsequently, I will first briefly discuss the case of a noninteracting many-fermion system 
(the Fermi gas), and then turn to Landau's theory of the interacting case (the liquid), from 
a phenomenological point of view and limited to the equilibrium properties. A much more 
detailed and complete exposition of these subjects can be found in the literature [|B|-[TC . 



A. The Fermi Gas 

In a noninteracting translationally invariant systems, the single-particle eigenstates are 
plane waves with energy eg = k 2 /2m. We will always use units so that h = 1. The ground 
state of an A-particle system is the well-known Fermi sea: all states up to the Fermi 
wavevector k F are filled, all the other states are empty. The energy of the last occupied state 
is the Fermi energy E F = k F 2 /{2m). 

It is usually convenient to define the Hamiltonian in a way so that the absolute ground 
state has a well-defined fixed particle number. This is achieved simply by including the 
chemical potential in the definition of the Hamiltonian, i.e. by writing 

# = (2-1) 

k 

where rig is the usual number operator, £g = e — /j,, and the spin summation is implicit. 
With this definition of the Hamiltonian, the elementary excitations of the Fermi gas are 

• addition of a particle at wavevector k (<5ng = 1). This requires \k\ > k F , and thus the 
energy of this excitation is eg — \i > 0. 

• destruction of a particle at wavevector k (<5ng = —1), i.e. creation of a hole. This 
requires \k\ < k F , and thus the energy is \i — eg > 0. 

These excitations change the total number of particles. Construction of states at constant 
particle number is of course straightforward: one takes one particle from some state k, with 
| A; | < k F , and puts it into a state k', with \k'\ > k F . These particle-hole excitations are 
parameterized by the two quantum numbers k,k' and thus form a continuum. The restric- 
tion on the allowed values of k, k' insures that all particle-hole states have positive energy. 
Higher excited states, i.e. states with many particles and many holes, are straightforwardly 
constructed, the only restriction being imposed by the Pauli principle. 

Thermodynamic quantities are easily obtained and are all determined by the density of 
states at the Fermi energy. For example, the specific heat obeys the well known linear law 
C{T) = 7 T, with 

7 = ^N(E F )k 2 B (2.2) 



3 



and similarly the (Pauli) spin susceptibility x an d the compressibility k are given by 



X = 2N(EM ,k = 2N(E f )/ P 2 . (2.3) 



It should be emphasized that eqs. Q2.2| ) to fl2.3| ) are valid for an arbitrary density of states, 
in particular in solids where bandstructure effects can change the electronic dispersion rela- 
tion quite drastically. Thus, for noninteracting electrons one expects the so-called "Wilson 
ratio" Rw = tt 2 fc|x/(3//|7) to be unity, independently of details of the bandstructure. Any 
deviation from unity is necessarily an indication of some form of interaction effect. 



B. Landau's theory: phenomenological 

1. Basic hypothesis 

Landau's theory is to a large extent based on the idea of a continuous and one-to-one 
correspondence between the eigenstates (ground state and excited states) of the noninter- 
acting and the interacting system. For this to be an acceptable hypothesis it is crucial that 
the interactions do not lead to any form of phase transition or symmetry-broken ground 
state. The standard way to represent this situation is the so-called "adiabatic switching 
on" of an interaction, parameterized by a variable V. On can easily imagine to numerically 
calculate eigenvalues of a Hamiltonian as a function of V (though this may be very difficult 
in practice). The basic hypothesis of Landau's theory then is that the ground state energy 
as well as the energy of the low-lying states as a function of V do not have any singularities 
due to some form of symmetry breaking. Moreover, it should be possible to label the low- 
lying states by the same quantum numbers as those of the noninteracting system, i.e. the 
number of excited particles and holes and their respective momenta. 

In particular one can consider a state obtained by adding a particle (with momentum 
\p\ > k-p) to the noninteracting ground state: 

\p,N + l) = a±\0,N) . (2.4) 

Here at is a fermion creation operator for momentum state p, and |0, N) is the iV-particle 
ground state of the noninteracting system. Now we add some form of particle-particle 
interaction. In a translationally invariant system, interactions conserve total momentum, 
and thus even after switching on the interaction the state still has total momentum p. 
However, the interaction of the added particle with the filled Fermi sea, as well as the 
interaction of the particles in the sea amongst themselves, will change the distribution of 
particles in /c-space, and of course also modify the energy of our state. The complex formed 
by the particle added at p and the perturbed distribution of the other particles is called a 
Landau quasiparticle. The Pauli principle implied \p\ > hp in the absence of interactions, 
and by the continuity hypothesis the same restriction remains valid in the interacting case. 
In particular, the value of hp, which imposes a lower limit on the allowed momentum of the 
quasiparticle, is unchanged by the interactions. 

Analogous considerations can be performed for a state obtained by destruction of a 
particle (e.g. creation of a hole): 

\p,N-l}=a_-\0,N) . (2.5) 
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Note that due to the momentum — p the total momentum of this state is indeed p. The state 
obtained after switching on the interactions (and which of course still obeys \p\ < kp) could 
be called a quasi-hole, but one generally refers to both quasi-particles and quasi-holes as 
quasi-particles. 

The quasi-particle concept has a certain number of limitations, mainly due to the fact 
that, as will be discussed below, the lifetime of a quasi-particle is finite. However, for 
excitations close to hp one has 1/r oc (e — E-p) 2 , i.e. the lifetime becomes much longer than 
the inverse excitation energy, and the quasi-particles therefore are reasonably well defined 
in the vicinity of the Fermi energy. In practice, this means that Landau's theory is useful for 
phenomena at energy scales much smaller than the Fermi energy, but inapplicable otherwise. 
In metals, where Ep rs 3 . . . 5eV, this restriction is not too serious when one is concerned 
with thermodynamic or transport properties. One should also note that the ground state 
energy itself has important contributions from states well below Ep, and therefore is not 
accessible to Landau's theory. On the other, the excitations above the ground state are the 
fundamental objects of Landau's theory. 



2. Equilibrium properties 

To be able to derive physical quantities from the picture of the low-energy excitations, 
one needs some information about the energetics of the quasiparticles and of their interac- 
tions. To be specific, starting from the ground state quasiparticle distribution 

n (k) = 1 if \k\ < kp 

= if \k\ > k F (2.6) 

one considers changes in quasiparticle occupation number of the form no(k) — > no(k)+Sn(k), 
i.e. Sn(k) = 1 represents an excited quasi-particle, Sn(k) = — 1 an excited quasi-hole (with 
the notation k = (k, a), and a =f , \ the spin index). The corresponding change in energy is 

SE = Y, £ °M k ) /(*> k')8n{k)8n{k') , (2.7) 

k ZU kk' 

where the first and second term represent the energy of a single quasi-particle and the 
interaction between quasiparticles, respectively. To be more precise, we assume that the 
chemical potential is included in the Hamiltonian, as in eq . ( | . 1|) . Consequently, vanishes 
on the Fermi surface, and, given that we are mainly interested in phenomena in the vicinity 
of kp, it is sufficient to retain the lowest order term in an expansion around \k\ = kp. One 
thus writes 

el = ^(\k\-kp) , (2.8) 
m 

thus defining the effective mass m* . The difference between the "bare mass m and m* is 
due to interaction effects which in principle can be calculated from a microscopic theory of 
the system in question. At the present level m* will be considered as a phenomenological 
parameter. 
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The energy of a quasi-particle added to the system is easily obtained from eq.( [2.7|) by 
calculating the difference in SE between a state with Sn(k) = 1 and a state with 5n(k) = 0. 
One finds 

e k = el + Wf^ k ') 6 < k ') » ( 2 - 9 ) 

i.e. the energy of an added quasi-particle is not just the "bare" quasiparticle energy eS 
but also depends, via the interaction term, on the presence or not of other quasi-particles. 
Given that the non-interacting particles obey Fermi-Dirac statistics, the quasi-particles do 
so too, and consequently, the occupation probability of a quasi-particle state is given by 

n(k) = - . (2.10) 

Note that the full quasi-particle energy and not the bare e° k enters this expression. In 
principle, n(k) thus has to be determined selfconsistently from eqs.( [2.9| ) and ( P-lOj ). 

For the subsequent calculations, it is convenient to transform the quasiparticle interaction 
f(k,k'). First, spin symmetric and antisymmetric /-functions are defined via 

fikik'D = f s {k,k') + nk,k') 

/(*T,*U) = f s (k,k')-f a (k,k') (2.11) 

Moreover, given the implicit restrictions of the theory, one is only interested in processes 
where all involved particles are very close to the Fermi surface. Under the assumption that 
the interaction functions are slowly varying as a function of k, one then can set \k\ = \k'\ = 
kp. Because of rotational symmetry, the /-functions then can only depend on the angle 
between k and k', called 9. One can then expand the /-function in a Legendre series as 

-* -* 00 k ■ P 

r s (k,k') = J2f a L S PL(cos6) , COS 6 = —^ , (2.12) 

L=0 fc F 

where the Pl are the Legendre polynomials. This formula can be inverted to give 
2 T A- 1 r 

fi S = —r^ / d 2 fiP L (cos#)r' s (M') • (2-13) 



47T 

Finally, one usually puts these coefficients into dimensionless form by introducing 

p a,s = (2 M) 

We are now in a position to calculate some equilibrium properties. The first one will be 
the specific heat at constant volume 

c " - hw (2 - 15) 

where U is the internal energy. The temperature-dependent part of U comes from thermally 
excited quasi-particles, as determined by the distribution ( |2.10| ). In principle, in this expres- 
sion Sk is itself temperature-dependent, because of the temperature dependent second term 
in eq.( p.9|) . However, one can easily see that this term only gives contributions of order T 2 , 
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and therefore can be neglected in the low-temperature limit. Consequently, one can indeed 
replace by e% and then one only has to replace the bare mass by m* in the result for a 
non-interacting system to obtain 



Cr 



m*k 



F;„2 



kiT 



(2.16) 



The spin susceptibility (at T = 0) is related to the second derivative of the ground state 
energy with respect to the (spin) magnetization M: 



X 



n 



d 2 E 
DM 2 



(2-17) 



At finite temperature the ground state energy Eq has to be replaced by the free energy 
F(T,N). Spin magnetization is created by increasing the number of f spin particles and 
decreasing the number of j spins (M = /j,b(N^ — N±)), i.e. by changing the Fermi wavevectors 
for up and down spins: kp — ► kp + 5kp for a =] and hp —>■ kp — 5kp for a =|. One has 



Sn(k, T) 

HK I) 

and the magnetization is 



1 if k F < \k\ < k F + 5k F 
— 1 if fcp — 8kp < \k\ < kp 



M 



n^skp 



7T' 



(2.18) 



(2-19) 



Using the 5n's from eq.( p.!8| ), the contributions from the first and second term in eq. ( p.7| ) 
are straightforwardly calculated as 



5 Ek = 



k F 3 
2TT 2 m* 



5kp 2 



5Eo = Q 



k F 3 



F"5kp 2 



(2.20) 



The coefficient Fq appears because the distortion Q2.18D is antisymmetric in the spin index 
and has spherical spatial symmetry (L = 0). From eqs.( |2.17| ) and ( p.20| ) one obtains the 
susceptibility as 



X 



jj^kpm* 



l + F$ 



7T- 



(2-21) 



Note that here (and contrary to the specific heat) interactions enter not only via m* but 
also explicitly via the coefficient Fq. The Wilson ratio then is 



R 



1 + FZ 



[2.22) 



Beyond the concrete result obtained here, it should be noted that the contributions from 
the two terms in eq.( |2.7| ) are both of order Skp 2 , even though formally they seem to be of 
order 8n and 8n 2 , respectively. This is in fact due to the "extra smallness" coming from the 
vanishing of e®* at the Fermi surface, whereas f(k, k') remains of course finite. On the other 
hand, it should also be clear that adding terms of order Sn 3 to eq. ( |2.7| ) would only lead to 
higher order corrections. 
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Following a similar reasoning, one can obtain a result for the compressibility of a Fermi 
liquid. The general definition of k is 



K = — — 



1 dQ 

nop 



d 2 E 

on 2 



[2.23) 



where at finite temperature one again has to replace the ground state energy by the free 
energy. For the calculation it is convenient to replace the variation of the volume at constant 
particle number by a corresponding variation of the particle number at constant volume. 
One then obtains an equivalent expression: 



2 d 2 e 
' dp 2 



(2.24) 



where p is the particle density and eo is the ground state energy density. The calculation 
then is very similar to the one for the spin susceptibility, only now hp is changed in the same 
way for up and down spins, and consequently F£ appear in the final result 



7rV(l 



(2.25) 



It is also interesting that in a translationally invariant system as we have considered 
here, the effective mass is not independent of the interaction coefficients. This can be seen 
the following way. Because of Galilean invariance a system moving at constant velocity v 
has an extra kinetic energy Nmv 2 /2. On the other hand, such a system can be represented 
by taking some particles from the left side of the Fermi sea to the right side, i.e. by a 5n(k) 
of dipolar symmetry. One can then use this Sn(k) to calculate the corresponding change in 
energy and compare to Nmv 2 /2. From the comparison one finds 



m 



l + F{/3 . 



[2.26) 



In the present phenomenological framework, the ^-coefficients are phenomenological 
parameters, to be determined from experiment. The standard example is liquid helium 
3, for which one finds [|l|,0 m*/m « 3, F S « 9, F{ « 5, F a « -0.7, Ff « -0.55. 
These values indicate rather strong interaction effects. Note the negative value of F£, which 
corresponds to a strong enhancement of the spin susceptibility. 



C. Conclusion 



The Landau theory of equilibrium properties is of rather limited quantitative predictive 
power because there is no prediction about the actual values of the Landau parameters. Its 
principal importance is of conceptual nature: from the single hypothesis about the existence 
of quasiparticles it follows that the low temperature equilibrium properties are very similar to 
those of a noninteracting system. In particular, the temperature dependencies are unaffected, 
only prefactors are interaction-dependent. Actual quantitative predictions are obtained 
when one extends the theory to nonequilibrium properties, using a Boltzmann equation. 
A new phenomenon predicted (and actually observed in 3 He fll3f ) is the existence of 
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collective excitations, called "zero sound" . This approach also allows the calculation of the 
quasiparticle lifetime and its divergence as the Fermi energy is approached, as well as the 
treatment of a number of transport phenomena. 

It does not seem generally possible to derive Landau's theory starting from some mi- 
croscopic Hamiltonian, apart possibly in perturbation theory for small interactions. It is 
however possible to formulate the basic hypotheses in terms of microscopic quantities. In 
particular, in the microscopic language the existence of quasiparticles translates into the 
existence of the quasiparticle pole in the one particle Green function: 

Q(k,iu n ) = - Zk . (2.27) 

Here zj: is the so-called quasiparticle weight, and gives rise to a jump in the momentum 
distribution function at kp of height zz, rather than unity in the noninteracting case. Simi- 
larly, the Landau interaction parameters can be related to two-particle vertex parts in the 
limit of vanishing momentum transfer. It should however be emphasized that this line of 
reasoning provides a microscopic interpretation of Landau's picture, rather than proving its 
correctness. Similar remarks apply to the calculated diverging quasiparticle lifetime: this at 
best shows that Landau's picture is internally consistent. 

As already mentioned, the ideal system for the application of Landau's theory is 3 He, 
which has both short-range interaction and is isotropic. The application to electrons in 
metals is more problematic. First, the interactions are long-ranged (Coulombic). This can 
however be accommodated by properly including screening effects. More difficulties, at least 
at the quantitative level, arise because metals are naturally anisotropic. This problem is not 
of fundamental nature: even when the Fermi surface is highly anisotropic, an expansion like 
eq. (|2.7|) can still be written down and thus interaction parameters can be defined. However, 
a simple Legendre expansion like eq.(|2.12|) is not in general possible, i.e. the description of 



the quasiparticle interaction in terms of a few parameters becomes impossible. An exception 
cases with a very nearly spherical Fermi surface like the alkali metals, where a determination 
of Landau parameters can indeed be attempted. It should be noticed that the difficulties 
with the Landau description of metals are not of conceptual nature and in particular do 
not invalidate the quasiparticle concept but are rather limitations on the usefulness of the 
theory for quantitative purposes. 



III. RENORMALIZATION GROUP FOR INTERACTING FERMIONS 

In this chapter, we will consider properties of interacting fermions in the framework 
of renormalization group theory. This will serve two purposes: first, the treatment one- 
dimensional interacting fermions, which will be considered in detail in the following chapters, 
gives rise to divergences which can only be handled by this approach. Results obtained in 
this way will be an essential ingredient in the subsequent discussion of "Luttinger liquids" . 
More generally, the renormalization group method will clarify the status of both Landau's 
Fermi liquid theory and the Luttinger liquid picture as renormalization group fixed points, 
thus establishing a link with a number of other phenomena in condensed matter physics. We 
will formulate the problem in terms of fermion functional integrals, as done by Bourbonnais 
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in the one-dimensional case |14| and more recently for two and three dimensions by Shankar 



15|J16|| . For the most part, I will closely follow Shankar's notation. 

Before considering the interacting fermion problem in detail, let us briefly recall the 
general idea behind the renormalization group, as formulated by Kadanoff and Wilson: one 
is interested in the statistical mechanics of a system described by some Hamiltonian H. 
Equilibrium properties then are determined by the partition function 

Z= £ e"^ = £ e"* , (3.1) 
configurations configurations 

where the second equality defines the action S = f3H. Typically, the action contains degrees 
of freedom at wavevectors up to some cutoff A, which is of the order of the dimensions of 
the Brillouin zone. One wishes to obtain an "effective action" containing only the physically 
most interesting degrees of freedom. In standard phase transition problems this is the 
vicinity of the point k = 0, however, for the fermion problem at hand the surface \k\ = kp is 
relevant, and the cutoff has to be defined with respect to this surface. In order to achieve 
this one proceeds as follows: 

1. Starting from a cutoff-dependent action 5(A) one eliminates all degrees of freedom 
between A and A/s, where s is a factor larger than unity. This gives rise to a new 
action S'(A' = A/s). 

2. One performs a "scale change" k — > sk. This brings the cutoff back to its original value, 
i.e. one obtains a new action S'(A). Because of the degrees of freedom integrated out, 
coupling constants (or functions) are changed. 

3. One chooses a value of s infinitesimally close to unity: s — 1 +e, and performs the first 
two steps iteratively. This then gives rise to differential equations for the couplings, 
which (in favorable circumstances) can be integrated until all non-interesting degrees 
of freedom have been eliminated. 



A. One dimension 

The one-dimensional case, which has interesting physical applications, will here be 
mainly used to clarify the procedure. Let us first consider a noninteracting problem, e.g. a 
one-dimensional tight-binding model defined by 

H = J2& a l a k > f* = -Ztcosk-n , (3.2) 

k 

where t is the nearest -neighbor hopping integral. We will consider the metallic case, e.g. the 
chemical potential is somewhere in the middle of the band. Concentrating on low-energy 
properties, only states close to the "Fermi points" ±fcp are important, and one can then 
linearize the dispersion relation to obtain 

H — VF^rk — k F )al r a kr , (3.3) 

fe,r=± 
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where vp = 2tsmkp is the Fermi velocity, and the index r differentiates between right- and 
left-going particles, e.g. particles close to k F and — k F . To simplify subsequent notation, we 
(i) choose energy units so that vp = 1, (h) translate fc-space so that zero energy is at k = 0, 
and (iii) replace the fc-sum by an integral. Then 

/A Ah 
A -rkal(k)a r (k) . (3.4) 

For the subsequent renormalization group treatment we have to use a functional integral 
formulation of the problem in terms of Grassmann variables (a detailed explanation of this 



formalism is given by Negele and Orland fll7|). The partition function becomes 



Z(A) = J X?0e~ 5(A) , (3.5) 

where U<p indicates functional integration over a set of Grassmann variables. The action is 

rP ( r A dk 1 
S(A) = J Q drlj2J_ A — 4>t{k,T)d T( j> r {k,T)+H{4>\4>)\ , (3.6) 

where the zero-temperature limit (3 — > oo has to be taken, and if(0*,0) indicates the 
Hamiltonian, with each a) replaced by a 0*, and each a replaced by a 0. Fourier transforming 
with respect to the imaginary time variable 

<f> r {k,T) =Tj2Mk^n)e-' luJnT K = 27r(n + 1/2)T) (3.7) 
and passing to the limit T — > one obtains the noninteracting action 

/OO f\(jj /"A Air 
— / —fi(k,u>)[-iLj + rk]<l> r (k,u,) . (3.8) 
-oo 271 J-A iTX 

We notice that this is diagonal in k and u which will greatly simplify the subsequent treat- 
ment. Because of the units chosen, uj has units of (length) -1 (which we will abbreviate as 
L _1 ), and then <fi r (k,uj) has units L 3 / 2 . 

We now integrate out degrees of freedom. More precisely, we will integrate over the strip 
A/s < \k\ < A, — oo < uj < oo. The integration over all uj keeps the action local in time. 
One then has 

Z{A) = Z(A,A/s)Z(A/s) , (3.9) 

where Z(A, A/s) contains the contributions from the integrated degrees of freedom, and 
Z(A/s), So (A/s) are of the form of eqs.( [3.5|) and (|3.8|) . The diagonality of So (A) leads to 
the factorized form of eq. (|3.9|) . Introducing the scale change 

k' = ks , J = us , 0' = 0s" 3/2 (3.10) 

one easily finds that in fact So (A/s) = So (A), i.e. the action does not change under scale 
change (or renormalization) , we are at a fixed point. One should notice that the scale change 
of k implies that k! is quantized in units of Ak' = 2ns /L, i.e. eliminating degrees of freedom 
actually implies that we are considering a shorter system, with correspondingly less degrees 
of freedom. This means that even though the action is unchanged the new Z(A) is the 
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partition function of a shorter system. To derive this in detail, one has to take into account 
the change in the functional integration measure intervening due to the scale change on 0. 

Before turning to the problem of interactions, it is instructive to consider a quadratic 
but diagonal perturbation of the form 

SS* = £ r T- t ^tik,u)ti.{k,u>)Mk,t»>) ■ (3-11) 

r= ±J-oc Z7T J -A 271 

We assume that fi(k,u) can be expanded in a power series 

fi(k, u) = /i 00 + fiiok + fx 01 iuj + ... (3. 12) 

Under the scale change ( |3 .10 ) one then has 

f^nrn * S f^nm • (3.13) 

There now are three cases: 

1. a parameter m nm grows with increasing s. Such a parameter is called relevant. This 
is the case for /z o- 

2. a parameter remain unchanged (/^io^oi)- Such a parameter is marginal. 

3. Finally, all other parameter decrease with increasing s. These are called irrelevant. 

Generally, one expects relevant parameter, which grow after elimination of high-energy 
degrees of freedom, to strongly modify the physics of the model. In the present case, the 
relevant parameter is simply a change in chemical potential and this doesn't change the 
physics much (the same is true for the marginal parameters). One can easily see that 
another relevant perturbation is a term coupling right- and left-going particles of the form 
m(0^0 2 + 020i)- This term in fact does lead to a basic change: it leads to the appearance 
of a gap at the Fermi level. 

Let us now introduce fermion-fermion interactions. The general form of the interaction 
term in the action is 

Sj= [ «(1234)0*(1)0*(2)0(3)0(4) . (3.14) 

J kcu 

Here 0(3) is an abbreviation for (f) r3 (ks, Us) , and similarly for the other factors, u is an 
interaction function to be specified. The integration measure is 



ik 



00 da;,- r A dk 



ku \JJ 2tt J-A 27T / 

xS(ki + fc 2 — k 3 — fc 4 )5(a;i + u 2 — a; 3 — a; 4 ) ... . (3.15) 

We now note that the dimension of the integration measure is L~ 6 , and the dimension of 
the product of fields is L 6 . This in particular means that if we perform a series expansion 
of u in analogy to eq.( 3.12|) the constant term will be s-independent, i.e. marginal, and all 
other terms are irrelevant. In the following we will thus only consider the case of a constant 
(k- and ^-independent) u. 
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These considerations are actually only the first step in the analysis: in fact it is quite 
clear that (unlike in the noninteracting case above) integrating out degrees of freedom will 
not in general leave the remaining action invariant. To investigate this effect, we use a more 
precise form of the interaction term: 

Si = 9i I EC + (l)^(2)^ + (3)t(4) 

J kui i 
ss' 

+92 f £0: + (l)0:^2)0 s ,_(3)0 s+ (4) . (3.16) 

J kuj i 
ss' 

Here we have reintroduced spin, and the two coupling constants g% and g 2 denote, in the 
original language of eq. (|3~3|) , backward ((kp, —kp) — > (— kp, kp)) and forward ((kp, —kp) — > 
(kp, —kp)) scattering. Note that in the absence of spin the two processes are actually iden- 
tical. 

Now, the Kadanoff-Wilson type mode elimination can be performed via 

e" 5 ' = J V(f)e- S , (3.17) 

where Dip denotes integration only over degrees of freedom in the strip A/s < \k\ < A. 
Dividing the field (f> into <f) (to be eliminated) and <$' (to be kept), one easily sees that the 
noninteracting action can be written as So = Sq((/)') +So(4>). For the interaction part, things 
are a bit more involved: writing 



Si — 22 Si t i — Si t o + Si , (3.18) 

where Sj t i contains i factors 4> we obtain 

e -S' = e -s»W-5/,o [ V4>e~ So ^~ Sl . (3.19) 



Because Sj contains up to four factors 0, the integration is not straightforward, and has to 
be done via a perturbative expansion, giving 



J V(j)e- So ^- §I = Z (A,A/s)exp 



oo ^ 

5^ ~}(Sj)o,con 



(3.20) 



where the notation (. . .)o, CO n indicates averaging over and only the connected diagrams 
are to be counted. The first order cumulants only give corrections to the energy and the 
chemical potential and are thus of minor importance. The important contributions come 
from the second order term (S] 2 )o,con which after averaging leads to terms of the form 
(p'* (p'* <fi' '(f)' , i.e. to corrections of the interaction constants g± t 2- The calculation is best done 
diagrammatically, and the four intervening diagram are shown in fig.|l|. One can easily see 
that not all of these diagrams contribute corrections to g\ or g 2 . Specifically, one has 

Sgi oc gig 2 [(a) + (c)] + 2gj(d) 

$92 oc (gl + g 2 2 )(a)+gl(b) (3.21) 
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FIG. 1. The diagrams intervening m the renormahzation 01 the coupling constants g\ and 
Note that in (b) the direction of one arrow is reversed with respect to (a), i.e. this is a particle-hole 
diagram. 



where the factor 2 for diagram (d) comes from the spin summation over the closed loop. 
Because the only marginal term is the constant in w(1234), one can set all external energies 
and momenta to zero. The integration over the internal lines in diagram (a) then gives 

dk ,<L, 1 1 ,a d*l 1 



In J 2n iu) — k —iuo — k J a/s 2n k 2n 

where s = 1 + d£, and similarly the particle-hole diagrams (b) to (d) give a contribution 
—d£/(2n). Performing this procedure recursively, using at each step the renormalized cou- 
plings of the previous step, one obtains the renormahzation group equations 

d9i 1 2 //* dgi 1 2 

- d J = -n 9 ^ ) ' -dl = -^ {e) ' (3 - 23) 

where s = e £ . These equations describe the effective coupling constants to be used after 
degrees of freedom between A and Ae e have been integrated out. As initial conditions one 
of course uses the bare coupling constants appearing in eq. fl3.16Q . 

Equations ( |3.23| ) are easily solved. The combination gi — 2g 2 is ^-independent, and one 
has further 

M = txS • ( 3 - 24 ) 

There then are two cases: 

1. Initially, g\ > 0. One then renormalizes to the fixed line g\ = gi(£ — > oo = 0, 
g\ — g2 — 9x1 i- e - one °f t ne couplings has actually vanished from the problem, but 
there is still the free parameter g\. A case like this, where perturbative corrections 
lead to irrelevancy, is called "marginally irrelevant" . 

2. Initially, g\ < 0. Then gi diverges at some finite value of I. One should however notice 
that well before the divergence one will have left the weak-coupling regime where 
the perturbative calculation leading to the eq. (|3.23| ) is valid. One should thus not 



overinterpret the divergence and just remember the renormalization towards strong 
coupling. This type of behavior is called "marginally relevant" . 

We will discuss the physics of both cases in detail in the next section. 

Two remarks are in order here: first, had we done a straightforward order-by-order 
perturbative calculation, integrals like eq. (|3.22| ) would have been logarithmically divergent, 
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both for particle-particle and particle-hole diagrams. This would have lead to inextricably 
complicated problem already at the next order. Secondly, for a spinless problem, the factor 
2 in the equation for gi(£) is replaced by unity. Moreover, in this case only the combination 
gi — gi is physically meaningful. This combination then remains unrenormalized. 



B. Two and three dimensions 



We will now follow a similar logic as above to consider two and more dimensions. Most 
arguments will be made for the two-dimensional case, but the generalization to three di- 
mensions is straightforward. The argument is again perturbative, and we thus start with 
free fermions with energy 



K 2 

£g = — - fi = v F k + 0(k 2 ) (v F = kp/m) . 



(3.25) 



We use upper case momenta K to denote momenta measured from zero, and lower case 

— # 

to denote momenta measured from the Fermi surface: k = \K\ — kp. The Fermi surface 
geometry now is that of a circle as shown in fig.^j. One notices in particular that states are 







' k F 











FIG. 2. Fermi surface geometry in two dimensions. 

now labeled by two quantum numbers which one can take as radial (k) and angular (8). 
Note that the cutoff is applied around the low-energy excitations at \K\ 
K = 0. The noninteracting action then takes the form 

So = kp r %l r p. f A ^(ke^i-iu - k)<p{kduj)] . 

J-oo Z7T JO Z7T J-A Z7T 



fcp, not around 



(3.26) 



One notices that this is just a (continuous) collection of one-dimensional action functional, 
parameterized by the variable 9. The prefactor kp comes from the two-dimensional integra- 
tion measure d 2 K = (kp + k)dkd6, where the extra factor k has been neglected because it 
is irrelevant, as discussed in the previous section. 

The general form of the interaction term is the same as in the one-dimensional case 



Si 



Kuj 



«(1234)0*(1)0*(2)0(3)0(4) 



(3.27) 
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however, the integration measure is quite different because of two-dimensional i?-space. 
Performing the integration over K4 and in the two-dimensional analogue of eq.( |3.15| ), 
the measure becomes 



Kuj 



M ^ f ^ r d|^^ e(A _ N) ... (3 . 28) 



\27ry y^J-oo 2tt Jo 2tt J -a 2tt 

Here K4 = K 1 + K 2 — K 3 . Now the step function poses a problem because one easily 
convinces oneself that even when K 12t3 are on the Fermi surface, in general K 4 can be far 
away from it. This is quite different from the one-dimensional case, where everything could 
be (after a trivial transformation) brought back into the vicinity of k — 0. 

To see the implications of this point, it is convenient to replace the sharp cutoff in 
eq.( |3.28| ) by a soft cutoff, imposed by an exponential: 

9(A- \k 4 \) ->exp(-|fc 4 |/A) • (3.29) 

Introducing now unit vectors f2j in the direction of Aj via Aj = (kp + ki)fli one obtains 

k 4 = \k F (&i + &2~ ^3) + htii + . . . I - k F « k F (\K\ - 1) 

A = fix +Q 2 - n 3 . (3.30) 

Now, integrating out variables leaves us with A — > A/s in eq.( p.2^ ) everywhere, including 
the exponential cutoff factor for /c 4 . After the scale change ( |3.10|) the same form of the 
action as before is recovered, with 

u'iK, u/i, 6[) = e^-^^-Mk/is, uJi/s, 0i) . (3.31) 

We notice first that nothing has happened to the angular variable, as expected because it 
parameterizes the Fermi surface which is not affected. Secondly, as in the one-dimensional 
case, the k and u> dependence of u is scaled out, i.e. only the values u(0, 0, 9i) on the 
Fermi surface are of potential interest (i.e. marginal). Thirdly, the exponential prefactor 
in eq.( 3.31| ) suppresses couplings for which |A| 7^ 1. This is the most important difference 



with the one-dimensional case. 

A first type of solution to |A| = 1 is 

Qi = ^3 => Q2 = ^4 , or 

Q ± = Q 4 =>■ Q 2 = Q 3 . (3.32) 

These two cases only differ by an exchange of the two outgoing particles, and consequently 
there is a minus sign in the respective matrix element. Both process depend only on the 
angle 812 between f2i and ^2, and we will write 

u(0, 0, 6 U 82, 0i, 8 2 ) = -u(0, 0, 8 1 , 82, 8 2 , 8,) = F{8 1 - 8 2 ) . (3.33) 

We can now consider the perturbative contributions to the renormalization of F. To lowest 
nontrivial (second) order the relevant diagrams are those of Fermi liquid theory and are 
reproduced in fig.|3]. Consider diagram (a). To obtain a contribution to the renormalization 
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(b) 

FIG. 3. Second order diagrams renormalizing 

the coupling function. Here P = Ki + K2,K3 = Ki — K, and K4 = Ki — K'. Q is the loop 
integration variable. 




FIG. 4. Phase space for diagram (a). The rings are the degrees of freedom to be integrates out 
between A and A/s. Note that only if P = are Q and P — Q simultaneously in the area to be 
integrates, giving a contribution of order d£. 



of F, both Q and P — Q have to lie in the annuli to be integrates out. As can be seen 
from fig.||, this will give a contribution of order d£ 2 and therefore does not contribute to a 



renormalization of F. The same is true (if we consider the first case in eq.( |3.32| )) for diagram 
(b). Finally, for diagram (c), because K is small, the poles of both intervening Green's 
functions are on the same side of the real axis, and here then the frequency integration gives 
a zero result. For the second process in eq. (|3.32 ) the same considerations apply, with the 
roles of diagrams (b) and (c) interchanged. The conclusion then is that F is not renormalized 
and remains marginal: 



dF 
~d£ 







(3.34) 



The third possibility is to have Qi 
Q% can be used to parameterize u: 



—Q2, ^3 = —^4- Then the angle between Q\ and 

u(o, 0, e 1 , -e u 3 , -63) = v(e 1 - e 3 ) . (3.35) 



In this case P = 0, and therefore in diagram (a) if Q is to be eliminated, so is —Q. Con- 
sequently, one has a contribution of order d£. For the other two diagrams, one finds again 
negligible contributions of order dl 2 . Thus, one obtains 



dV(9 1 - 9 3 



1 



8vr 2 Jo 2tt 



27r d9 

™v{(h - o)v(6 - e 3 ) 



(3.36) 
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(3.38) 



This is a renormalization equation for a function, rather than for a constant, i.e. one here 
has an example of a "functional renormalization group" . Nevertheless, a Fourier transform 

Vx = / TT^ e V{e) (3.37) 

JO 2,71 

brings this into a more standard form: 

dV x (t) _ Vx(i) 2 
d£ 4tt ' 

This has the straightforward solution 

From eqs.( |3.34D and eq.( |3.39| ) there are now two possibilities: 

1. At least one of the V\ is negative. Then one has a divergence of V\(£) at some fi- 
nite energy scale. Given that this equation only receives contributions from BCS-like 
particle-particle diagrams, the interpretation of this as a superconducting pairing in- 
stability is straightforward. The index A determines the relative angular momentum 
of the particles involved. 

2. All V\ > 0. Then on has the fixed point V\ = 0, F(6i — 82) arbitrary. What is the 
underlying physics of this fixed point? One notices that here one has 63 = 61,64 = 6 2 , 
i.e. the marginal term in the action is 4>g 1 4 > g 2 4>g 1 4>0 2 - In the operator language, this 
translates into 



l int 



(Wid^n^nft, . (3.40) 



We now can recognize this as an operator version of Landau's energy functional, 
eq.(2.7|). The fixed point theory is thus identified as Landau's Fermi liquid theory. 



One can notice that the fixed point Hamiltonian, eq. ( |3.4U , has a very large symmetry group: 
given that n k = a\a k , this term (as well as the kinetic energy term that only contains n,&) is 
invariant under the U(l) transformation 

a k -> e*>a k (3.41) 

with arbitrarily /c-dependent ipk, i.e. the symmetry group is XJ(1) N , with N the number of 
points on the Fermi surface. On the other hand, a standard interaction Hamiltonian has the 
form 

Him = ^qY, V (v) a l+ q a k>-q a k> a k , ( 3 - 42 ) 
kk'q 

and this is only invariant under a global U(l) transformation. 

The generalization of the above to three dimensions is rather straightforward. In addition 
to the forward scattering amplitudes F, scattering where there is an angle 0i2 ; 34 spanned by 
the planes (f2i, f2 2 ) and (0 3 , O4) is also marginal. For 0i2 ; 34 7^ these processes give rise to 
a finite quasiparticle lifetime, however they do not affect equilibrium properties. The (zero 
temperature) fixed point properties thus still only depend on amplitudes for 0i2 ; 34 = 0, i.e. 
the Landau /-function. 
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IV. LUTTINGER LIQUIDS 



The Fermi liquid picture described in the preceding two sections is believed to be rel- 
evant for most three-dimensional itinerant electron systems, ranging from simple metals 
like sodium to heavy-electron materials. The most studied example of non-Fermi liquid 
properties is that of interacting fermions in one dimension. This subject will be discussed 



in the remainder of these lecture notes. We have already started this discussion in sec. Ill A 



where we used a perturbative renormalization group to find the existence of one marginal 
coupling, the combination g 1 — 2g 2 . This approach, pioneered by Solyom and collaborators 
in the early 70's, can be extended to stronger coupling by going to second or even third 
order [|l||] in perturbation theory. A principal limitation remains however the reliance on 
perturbation theory, which excludes the treatment of strong-coupling problems. An alterna- 
tive method, which allows one to handle, to a certain extent, strong-interaction problems as 
well, is provided by the bosonization approach, which will be discussed now and which form 
the basis of the so-called Luttinger liquid description. It should be pointed out, however, 
that frequently entirely analogous results can be obtained by many-body techniques. |20|-|22 



A. Bosonization for spinless electrons 



The bosonization procedure can be formulated precisely, in the form of operator identi- 
ties, for fermions with a linear energy-momentum relation, as discussed in section |111 A| . To 
clarify notation, we will use a-(6-)operators for right- (left-) moving fermions. The linearized 
noninteracting Hamiltonian, eq. (|3.3| ) then becomes 



H = v F ^{(k ~~ k F )a\a k + 

k 



-k - k F )b[b k } 



(4.1) 



and the density of states is N(Ep) = 1/ (kvf)- In the Luttinger model, p3|,f4|| one generalizes 
this kinetic energy by letting the momentum cutoff A tend to infinity. There then are two 
branches of particles, "right movers" and "left movers" , both with unconstrained momentum 
and energy, as shown in fig.|5]. At least for weak interaction, this addition of extra states 




FIG. 5. Single-particle energy spectrum of the Luttinger model. Occupied states are shown in 
grey, the dark grey area represents the states added to make the model solvable. 
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far from the Fermi energy is not expected to change the physics much. However, this 
modification makes the model exactly solvable even in the presence of nontrivial and possibly 
strong interactions. Moreover, and most importantly, many of the features of this model 
carry over even to strongly interacting fermions on a lattice. 

We now introduce the Fourier components of the particle density operator for right and 
left movers: 



p+(q) = Y,' 



l k+q 



a k 



(4.2) 



The noninteracting Hamiltonian (and a more general model including interactions, see be- 
low) can be written in terms of these operators in a rather simple form and then be solved 
exactly. This is based on the following facts: 

1. the density fluctuation operators p± obey Bose type commutation relations: 



[p+(-q)>p+(q')\ = \p-(q),p-(-q')]= l 

\p + (q), P 4q')} = o . 



qh 
qq '2n 



(4.3) 
(4.4) 



The relation Q4.4Q as well as eq.( [4.3| ) for q ^ q' can be derived by straightforward 
operator algebra. The slightly delicate part is eq . ([13] ) for q = q'. One easily finds 



[p+(-q),p+(q)} 



k 



(4.5) 



where h k is an occupation number operator. In a usual system with a finite interval of 
states between — fcp and kp occupied, the summation index of one of the n operators 
could be shifted, giving a zero answer in eq. (|4.5|) . In the present situation, with an 
infinity of states occupied below kp (see fig.|6|), this is not so. Consider for example the 
ground state and q > 0. Then each term in eq.(|4.5|) with k F < k < k F + q contributes 



(a) 



(b) 



(O 







F u 

FIG. 6. Occupied states (heavy line) for right moving fermions (a), left moving fermions (b), 
and in a standard model without extra states at negative energy (c). 



unity to the sum, all other terms vanish, thus establishing the result (|4.3| ). More 
generally, consider a state with all levels below a certain value ko (< k F ) occupied, but 
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an arbitrary number of particle hole pairs excited otherwise. One then has, assuming 
again q > 0, 



5^(n fc _, - h k ) = XI + H ~ "*) ( 4 - 6 ) 

A: yfc>fco k<ko J 

= H ~ "*) ( 4 - 7 ) 

k>ko 

= h k ~ h k (4.8) 

k>k() — q k>ko 

= £ «. = | ■ (4-9) 

feg— q<k<ko 

The result is independent of fco, and one thus can take the limit fco ~~ > — oo. Together 
with an entirely parallel argument for p_, this then proves eq.( [4.3|) . Moreover, for q > 
both p+(—q) and P-(q) annihilate the noninteracting groundstate. One can easily 
recover canonical Bose commutation relations by introducing normalized operators, 



e.g. for q > 6j = y2n/(qL)p + (q) would be a canonical creation operator, but we 
won't use this type of operators in the following. 

2. The noninteracting Hamiltonian obeys a simple commutation relation with the density 
operators. For example 

[H ,p+(q)]=v F qp + (q) , (4.10) 

i.e. states created by p+(q) are eigenstates of H , with energy vpq. Consequently, the 
kinetic part of the Hamiltonian can be re-written as a term bilinear in boson operators, 
i.e. quartic in fermion operators: 

#o = ^EMsKH?) + P-(-q)p-(q)] ■ (4.11) 

L q>0 

This equivalence may be made more apparent noting that p+(q) creates particle-hole 
pairs that all have total momentum q. Their energy is Sk+ q — £&, which, because of 
the linearity of the spectrum equals v F q, independently of k. Thus, states created by 
p+{q) are linear combinations of individual electron-hole excitations all with the same 
energy, and therefore are also eigenstates of (|4.1|) . 



3. The above point shows that the spectra of the bosonic and fermionic representations 
of H are the same. To show complete equivalence, one also has to show that the de- 
generacies of all the levels are identical. This can be achieved calculating the partition 
function in the two representations and demonstrating that they are equal. This then 
shows that the states created by repeated application of p± on the ground state form 
a complete set of basis states f24|.|2l 
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We now introduce interactions between the fermions. As long as only forward scattering 
of the type (k F ] —k F ) — > (k F ; —k F ) or (k F ; k F ) — > (k F ; k F ) is introduced, the model remains 
exactly solvable. The interaction Hamiltonian describing these processes takes the form 

H^t = ^^{2# 2 (g)p+(g)p-(-g) 

+9i(q)\p+{q)p+{-q)+P-(-q)P-(Q)]} ■ ( 4 - 12 ) 

Here, g 2 {q) and g±(q) are the Fourier transforms of a real space interaction potential, and in 
a realistic case one would of course have g 2 {q) = 9i(q), but it is useful to allow for differences 
between g 2 and g&. For Coulomb interactions one expects g 2 , g^ > 0. In principle, the long- 
range part of the Coulomb repulsion leads to a singular g-dependence. Such singularities in 
the gi can be handled rather straightforwardly and can lead to interesting physical effects 
as will be discussed below. Here I shall limit myself to nonsingular g 2l g±. Electron-phonon 
interactions can lead to effectively attractive interactions between electrons, and therefore 
in the following I will not make any restrictive assumptions about the sign of the constants. 
One should however notice that a proper treatment of the phonon dynamics and of the 
resulting retardation effects requires more care ||26| . 



Putting together Q4.ll ) and ( 4.12 ), the complete interacting Hamiltonian, the Tomonaga- 



Luttinger model, then becomes a bilinear form in boson operators that is easily diagonalized 
by a Bogolyubov transformation. A first consequence is the expression for the excitation 
spectrum 

"(ff) = \1\[{VF + <?4(g)/(27r)) 2 - (g 2 (q)/(2n)) 2 ] 1 / 2 . (4.13) 

The diagonal boson operators are linear combinations of the original p operators, and con- 
sequently, these elementary excitations are collective density oscillations, their energy being 
determined both by the kinetic energy term and the interactions. 

We note here that in order for the Bogolyubov transformation to be a well-defined 
unitary transformation, g 2 (q) has to decrease at large q at least as |g| -1 / 2 . on the other hand, 
the large-g behavior of g 2 is unimportant for the low-energy properties of the model. We 
therefore in the following will almost always use a g-independent g 2 and #4. An approximate 
and frequently used way to cure the divergences arising due to this procedure is to keep the 
parameter a in subsequent formulae as a finite short-distance cutoff, of the order of a lattice 
spacing. One can then also include the "backward scattering" (k F ;—k F ) — > (— k F ', k F ), 
because for spinless electron this is just the exchange analogue of forward scattering and 
does not constitute a new type of interaction. It is worthwhile emphasizing here that the 
solution is valid for arbitrarily strong interactions, no perturbative expansion is needed! 

Up to this point, the construction does not allow for a direct calculation of correlation 
functions like the one-particle Green function or more generally any function involving 
individual creation or destruction operators. This type of correlation function becomes 
tractable by representing single particle operators in terms of the boson operators. To this 
end, we introduce the field operators 

m = -7E -e- aW2 - ipx [pM + P-ip)\ - , (4.14) 
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U(x) 



1 



aM/2 - ipx [p + {p)-p^{p)] + J/L 
N = N,+N_ , J = N + -N_ . 



(4.15) 
(4.16) 



Here N± is the number of particles added to the ground state on the right- and left-moving 
branch. Because addition of a particle changes both N and J, one has the "selection rule" 
(— 1)^ = (— 1) J , and a is a cutoff parameter which (at least in principle, see the discussion 
above) has to be set to zero in the end of any calculation. and IT then obey canonical 
boson commutation relations: 



[0(e): !%)] = i6(x-y) , 
and 4> is related to the local particle density via 
dcj)/dx = -7r(p(x) - po) , 



(4.17) 



(4.18) 



where po is the average particle density in the ground state. More precisely, in a lattice 
model this would represent the slowly varying components (q ~ 0) of the density, whereas 
components with q « 2k F would correspond to crossproducts between ip±. 
The expression for the single fermion operators then is 



lim 



zU± exp [±ikpx =F i<fr(x) + i9{x)\ 



(4.19) 



where the upper and lower sign refer to right- and left-moving electrons respectively, and 

e{x) = 7t rn(x')dx' 



E 



-a\p\/2— ipx 



P 



\P+{P) -P-(P)] + J 



nx 

T 



(4.20) 



The [/-operators decrease the total particle number on one of the branches by unity and 
are necessary because the boson fields all conserve the total particle number. A detailed 
derivation of the important eq.( [4.19| ) as an operator identity is given in the literature p4] , |25 



However, a simple plausibility argument can be given: creating a particle at site x means 
introducing a kink of height 7r in 0, i.e. at points on the left of x </> has to be shifted by 
7i. Displacement operators are exponentials of momentum operators, and therefore a first 
guess would be ip(x) m exp(i7r U(x')dx'). However, this operator commutes with itself, 
instead of satisfying canonical anticommutation relations. Anticommutation is achieved by 
multiplying with an operator, acting at site x, that changes sign each time a particle passes 
through x. Such an operator is exp(±i0(x)). The product of these two factors then produces 

(Hi). 

The full Hamiltonian can also be simply expressed in terms of and n. Neglecting the 
momentum dependence of the one easily finds 



H = Hq + H, 



int 



dx 



TTUK , 



-n(x) 2 + 



u 



2ttK 



(4.21) 



This is obviously just the Hamiltonian of an elastic string, with the eigenmodes correspond- 
ing to the collective density fluctuations of the fermion liquid. It is important to notice that 
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these collective modes are the only (low-energy) excited states, and that in particular there 
are no well-defined single particle excitations, nor are there the incoherent particle-hole pair 
excitations typical of a Fermi gas. The parameters in ( 4.21| ) are given by 



u=[(v F + g 4 /(2n)) 2 -g 2 2 /(2n) 2 ] 1 / 2 , K 



2nv F + g A - g 2 



1/2 



27tv f + g 4 + g 2 _ 

The energies of the eigenstates are u(q) = u\q\, in agreement with eq. Q4.13Q . 



(4.22) 



From the continuity equation, the expression ( 4.18 ) for the local particle density and the 



equation of motion of the (number) current density is 

j(x) = uKU(x) . (4.23) 

Note in particular that for g 2 = g± one has uK = v F , i.e. the expression for the current 
density is interaction-independent. The relation uK = v F holds in particular for systems 
with full (Galilean) translational invariance. On the other hand, in the continuum limit of 
lattice systems this relation is in general not true. 

The most remarkable result here is the "collectivization" of the dynamics: there are no 
quasiparticle-like excitations. In fact there is a rather simple physical picture explaining this: 
imagine accelerating one particle a little bit in one direction. Very soon it will hit its neighbor 
and transmit its momentum to it, and the neighbor will in turn transmit its momentum to 
a further neighbor, and so on. Quite quickly, the initial localized motion will have spread 
coherently through the whole system. This picture can be formalized noting that in one 
dimension the difference between a solid and a fluid is not well-defined: whereas is higher 
dimensions solids and fluids are differentiated by the presence or absence of long-wavelength 
transverse modes, no transverse modes can exist in a system with movement along only one 
direction. The long-wavelength modes thus can equally well be considered as the phonons 
of a one-dimensional crystal. |27|,[2£| Note that on the contrary in dimensions larger than 
one the neighbors of any given particle can be pushed aside, giving rise to a backflow that 
allows the particle to move trough the system more or less freely. 

Rather than discussing the physics of the spinless case in detail, we will turn now to the 
more interesting case of fermions with spin. 

B. Spin 1/2 fermions 

In the case of spin-1/2 fermions, all the fermion operators acquire an additional spin 
index s. Following the same logic as above, the kinetic energy then takes the form 

Hq = v F J2{(k-k F )ala k>s + {-k-k F )blb k } 



2nv F 



E [p+A<i)p+A-<i) + p-,s(-i)p-M] > ( 4 -24) 



L q>Q,s 

where density operators for spin projections s =|, j have been introduced: 

P+M = J2 a l+ q , s ak,s , P-M = Y, b l+ q ,s b k,s ■ (4-25) 
k k 
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There are now two types of interaction. First, the "backward scattering" {k F , s; —k F ,t) 
— > {—k F ,s;k F ,t) which for s ^ t cannot be re- written as an effective forward scattering 
(contrary to the spinless case). The corresponding Hamiltonian is 



HintA — — ^2 9l a k,s^p,t a p+2k F +q,tb k _2k 
kpqst 



(4.26) 



And, of course, there is also the forward scattering, of a form similar to the spinless case 

1 ■ y £{292(q)p + AQ)P-A-Q) ( 4 -27) 



H 



int,2 



2L 



qst 



+9Aq)\p+Aq)p+A-q) + p-A~q)p-Aq)]} 



(4.28) 



To go to the bosonic description, one introduces <fr and II fields for the two spin projections 
separately, and then transforms to charge and spin bosons via PjCr = ± <fi±)/A2, H p (T = 
(rif ± IIj,)/Vz. The operators <p u and TI^ obey Bose-like commutation relations: 

[(l>u{.x),UAv)] = iSunS(x-y) , 

and single fermion operators can be written in a form analogous to ( |4.19| ): 

tp±A x ) 



lini —= 



exp 



±ik F x - i(±((j)p + s(f) a ) + {Op + sO a 



(4.29) 



where u {x) = 7r J x U u {x')dx'. 
The full Hamiltonian H = H Q 



H 



intA 



Hi n t 2 then takes the form 



H = H + H a + 



2<?i 



J dx cos{A8(f) cr ) . 



(2vr«) 2 

Here a is a short-distance cutoff, and for v = p, cr 



with 



K v 



dx 



tiu v K v 2 



it 



2ttK, 



Tr{d x 4>u) 



(4.30) 



(4.31) 



2tiv f + 2g4^ + g v 
2nv F + 2# 4 „ - g v 



-,1/2 



(4.32) 



and g p = gi—2g 2 , g a — g%, g±,p = g±, <74, CT = 0. For a noninteracting system one thus has u v = 
v F (charge and spin velocities equal!) and K v = 1. For g 1 = 0, ( |4.30| ) describes independent 
long-wavelength oscillations of the charge and spin density, with linear dispersion relation 
oj v {k) = u u \k\, and the system is conducting. As in the spinless case, there are no single- 
particle or single particle-hole pair excited states. This model (no backscattering) , usually 
called the Tomonaga-Luttinger model, is the one to which the bosonization method was 
originally applied p3l^,p9|. 
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For gi ^ the cosine term has to be treated perturbatively. We have already obtained the 
corresponding renormalization group equations in the previous section (eq.( |3.23 )). In par- 
ticular, for repulsive interactions (gi > 0), g% is renormalized to zero in the long-wavelength 
limit, and at the fixed point one has K* = 1. The three remaining parameters in ( f4.30|) then 
completely determine the long-distance and low-energy properties of the system. 

It should be emphasized that ( |4.30| ) has been derived here for fermions with linear energy- 
momentum relation. For more general (e.g. lattice) models, there are additional operators 
arising from band curvature and the absence of high-energy single-particle states. |24j] One 
can however show that all these effects are, at least for not very strong interaction, irrelevant 
in the renormalization group sense, i.e. they do not affect the low-energy physics. Thus, 
( U-3(\ ) is still the correct effective Hamiltonian for low-energy excitations. The lattice effects 
however intervene to give rise to "higher harmonics" in the expression for the single-fermion 
operators, i.e. there are low energy contributions at wavenumbers q ~ (2m + l)kp for 
arbitrary integer m. 



The Hamiltonian ( |4 . 3 0|) also provides an explanation for the physics of the case of neg- 
ative gi, where the renormalization group scales to strong coupling (eq.( |3.23 )). In fact, if 
\gi | is large in ( |4.30|) , it is quite clear that the elementary excitations of CT will be small 
oscillations around one of the minima of the cos term, or possibly soliton-like objects where 
</> CT goes from one of the minima to the other. Both types of excitations have a gap, i.e. for 
gi < one has a gap in the spin excitation spectrum, whereas the charge excitations remain 
massless. This can actually investigated in more detail in an exactly solvable case. |30 



1. Spin-charge separation 



One of the more spectacular consequences of the Hamiltonian (|4.30|) is the complete 



separation of the dynamics of the spin and charge degrees of freedom. For example, in 
general one has u a ^ u p , i.e. the charge and spin oscillations propagate with different 
velocities. Only in a noninteracting system or if some accidental degeneracy occurs does 
one have u a = u p = vf- To make the meaning of this fact more transparent, let us create 
an extra particle in the ground state, at t = and spatial coordinate xq. The charge and 
spin densities then are easily found, using p(x) = — (\/2/iT)d(f) p /dx (note that p(x) is the 
deviation of the density from its average value) and a z (x) = — (\/2/i{)d(i) a /dx : 

(Q\il) + {x )p(x)4) ] + {x Q )\0) = 8{x-xq) , 

(0\iJ) + (x )(r g (x)1>l(xa)\0) = 5{x - x ) . (4.33) 

Now, consider the time development of the charge and spin distributions. The time- 
dependence of the charge and spin density operators is easily obtained from fl4.3Up (using 
the fixed point value gi = 0), and one obtains 

(0\4> + (x )p(x,t)4>i(x )\0) = S(x-x -u p t) , 

(0|V + (x )a 2 (x,t)^(xo)|0) = 5(x - x - u a t) . (4.34) 

Because in general u a ^ u p , after some time charge and spin will be localized at completely 
different points in space, i.e. charge and spin have separated completely. An interpretation of 
this surprising phenomenon in terms of the Hubbard model will be given in sec.( [VIB"D . For 
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simplicity we have taken here K p = 1, otherwise there would also have been a left-moving 
part in the time-dependent density. 

Here a linear energy-momentum relation has been assumed for the electrons, and conse- 
quently the shape of the charge and spin distributions is time-independent. If the energy- 
momentum relation has some curvature (as is necessarily the case in lattice systems) the 
distributions will widen with time. However this widening is proportional to y/i, and there- 
fore much smaller than the distance between charge and spin. Thus, the qualitative picture 
of spin-charge separation is unchanged. 



2. Physical properties 



The simple form of the Hamiltonian ( }4.30|) at the fixed point g\ = makes the calculation 



of physical properties rather straightforward. The specific heat now is determined both by 
the charge and spin modes, and consequently the specific heat coefficient 7 is given by 

l/lo = ^( v F/u p + v F /u a ) . (4.35) 

Here 70 is the specific heat coefficient of noninteracting electrons of Fermi velocity Dp. 

The spin susceptibility and the compressibility are equally easy to obtain. Note that 
in ( |4.30|) the coefficient u a /K a determines the energy necessary to create a nonzero spin 
polarization, and, as in the spinless case, u p /K p fixes the energy needed to change the 
particle density. Given the fixed point value K* = 1, one finds 

x/Xo = v F /u a , k/kq = v f K p /u p , (4.36) 

where xo and k are the susceptibility and compressibility of the noninteracting case. 

The quantity n p (x) is proportional to the current density. As before, the Hamiltonian 
commutes with the total current, one thus has 

a(uj) = 2K p u p 5(u) + o reg (u) , (4.37) 

i.e. the product K p u p determines the weight of the dc peak in the conductivity. If the total 
current commutes with the Hamiltonian a reg vanishes, however more generally this part of 
the conductivity varies as u 3 at low frequencies. fl3l | 



The above properties, linear specific heat, finite spin susceptibility, and dc conductivity 
are those of an ordinary Fermi liquid, the coefficients u p ,u ai and K p determining renormal- 
izations with respect to noninteracting quantities. However, the present system is not a 
Fermi liquid. This is in fact already obvious from the preceding discussion on charge-spin 
separation, and can be made more precise calculating the single-particle Green function 
using the representation fl4.29|) of fermion operators. One then obtains the momentum 
distribution function in the vicinity of hp 

n k ~ n k F — const. sign(A; — k F )\k — kp\ s , (4.38) 

and the single-particle density of states (i.e. the momentum-integrated spectral density): 

N(lo) « \co\ 5 . (4.39) 
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In both cases 5 = [K p + 1/K P — 2)/4. Note that for any K p ^ 1, i.e. for any nonvanishing 
interaction, the momentum distribution function and the density of states have power-law 
singularities at the Fermi level, with a vanishing single particle density of states at Ep. 
This behavior is obviously quite different from a standard Fermi liquid which would have 
a finite density of states and a step-like singularity in n^. The absence of a step at kp in 
the momentum distribution function implies the absence of a quasiparticle pole in the one- 
particle Green function. In fact, a direct calculation of the spectral function A(k,u) P2]J5B 



shows that the usual quasiparticle pole is replaced by a continuum, with a lower threshold 
at mxn(u u )(k — kp) and branch cut singularities at uj = u p p and uj = u a p: 

The coefficient K p also determines the long-distance decay of all other correlation func- 
tions of the system: Using the representation (|4.29|) the charge and spin density operators 
at 2kp are 

p2ikpx 

Ocdw(x) = = lim e- iv ^'cos[V2^] , (4.40) 

p 2ikpx 

OsdwM = E^-,^+,- s = lim e-^'cosfv^] • (4.41) 

Similar relations are also found for other operators. It is important to note here that all 
these operators decompose into a product of one factor depending on the charge variable 
only by another factor depending only on the spin field. Using the Hamiltonian ( 4.30|) at 



the fixed point g{ = one finds for example for the charge and spin correlation functions^ 

(n(z)n(O)) = K p /(txx) 2 + A l cos(2k F x)x- l - Kp \n-' i/2 (x) 

+ A 2 cos(Ak F x)x- iKp + ... , (4.42) 
(S(x)-S(O)) = 1/{tcx) 2 + B 1 cos{2k F x)x- 1 - Kp ln 1/2 {x) + . . . , (4.43) 



with model dependent constants A i} Bi. The ellipses in ( |4.42| ) and ( f4.43| ) indicate higher 



harmonics of cos(2kpx) which are present but decay faster than the terms shown here. 
Similarly, correlation functions for singlet (SS) and triplet (TS) superconducting pairing are 

(Ol s (x)O SS (0)) = Cx- 1 - 1 ^\n^ 2 (x) + ... , 
(4s>)CW0)> = Dx- 1 ~ 1 / K >\n 1 / 2 (x) + ... . (4.44) 



The logarithmic corrections in these functions |35| have been studied in detail recently 
P^|-pU[. The corresponding susceptibilities (i.e. the Fourier transforms of the above corre- 
lation functions) behave at low temperatures as 

Xcw(T)«rVi|l n ( T )|-^ ,Xsdw(T) ~ T^ -1 ! ln(T)! 1 / 2 , 

Xs S (T)™T 1 ' K >- 1 \H'n\-*' 3 ^Tsm^T'/^lHT)^ 2 , (4.45) 

i.e. for K p < 1 (spin or charge) density fluctuations at 2kp are enhanced and diverge at 
low temperatures, whereas for K p > 1 pairing fluctuations dominate. The "phase diagram", 
showing in which part of parameter space which type of correlation diverges for T — * is 
shown in fig.[7[ 



x The time- and temperature dependence is also easily obtained, see [34| 
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FIG. 7. Phase diagram for interacting spin-1/2 fermions. 



A remarkable fact in all the above results is that there is only one coefficient, K p , which 
determines all the asymptotic power laws, i.e. there are scaling relations between the ex- 
ponents of different correlation functions. These relations would have been impossible to 
obtain within a purely perturbative approach. [18] The correlation functions with their 



power law decay determine experimentally accessible quantities: the 2k f and Akp charge 
correlations lead to X-ray scattering intensities l2k F ~ T Kp , pa T AKp ~ l , and similarly 
the NMR relaxation rate due to 2kp spin fluctuations varies as 1/Tj ~ T Kp . Experimen- 



tal observations of the NMR relaxation rate j40] and of photoemission spectra j|T] on the 
quasi-one-dimensional organic compound (TMTSF) 2 PF 6 suggest K p pa 0.25, i.e. rather 
strong non-Fermi liquid behavior. It is currently unclear whether this interpretation is 
compatible with the transport properties. 

We here re-emphasize the two important properties of spin-1/2 interacting fermions in 
one dimension: (i) correlation functions show power-law decay, with interaction-dependent 
powers determined by one coefficient, K p ; and (ii) "spin-charge separation": spin and charge 
degrees of freedom propagate with different velocities. Both these properties invalidate the 
Landau quasiparticle concept in one dimension. Rather, the name "Luttinger liquid" has 
been coined to characterize the properties of one-dimensional interacting fermions. J24| 



3. Long-range interactions 

The above calculations can be straightforwardly generalized to the case of long-range in- 
teractions. [0 Of interest is the case of unscreened Coulomb interactions (V(r) = e 2 /r) 
which for example is of relevance for the physics of an isolated quantum wire. The 
short-distance singularity has to be cut off, and for example in a wire of diameter d 
an approximate from would be V(r) = e 2 /^/x 2 + d 2 , leading to a Fourier transform 
V(o) = 92(q) = 9i{.<i) = 2e 2 K (qd). The long-range nature of the interaction is only of 
importance for the forward-scattering processes, and these appear only in the charge part 
of the Hamiltonian which remains diagonal. The elementary excitations then are found to 
be charge oscillations (plasmons), with energy-momentum relation 

= + fifi)(l -9i + 2^(g))] V2 (4.46) 
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where V(q) = V(q)/ (ttvf) and g\ = gi/[2^VF)- The long-wavelength form, uj p (q) ~ 
|g 2 In gl 1 / 2 , agrees with RPA calculations |43|,|44]], however, the effect of g±, which is a short- 
range exchange contribution, is usually neglected in those calculations. The spin modes are 



still given by uj a (q) = u a \q\, with u a = vpJl — g\. Distinct charge- and spin modes have 



indeed been observed in Raman scattering from quantum wires. 

Correlation functions can be similarly calculated and one obtains for example 

(p(x)p(O)) = A 1 cos(2k F x) exp(— c 2 V\nx)/x 

+A 2 cos(4k F x)exp(-4c 2 Vhix) + ... , (4.47) 



where A lj2 are interaction dependent constants, c 2 = y(l + gi)irvF/e 2 , and only the most 
slowly decaying Fourier components are exhibited. The most interesting point here is the 
extremely slow decay (much slower than any power law!) of the 4k p component, showing 
an incipient charge density wave at wavevector Akp. This slow decay should be compared 
with the power-law decay found for short-range interactions (eq.( |4.42| )). The Akp oscillation 
period is exactly the average interparticle spacing, i.e. the structure is that expected for 
a one-dimensional Wigner crystal. Of course, because of the one-dimensional nature of 
the model, there is no true long-range order, however, the extremely slow decay of the 
Akp oscillation would produce strong quasi-Bragg peaks in a scattering experiment. It is 
worthwhile to point out that this Akp contribution arises even if the Coulomb interaction is 
extremely weak and depends only on the long-range character of the interaction. On the 
other hand, any 2kp scattering is considerably weaker, due to the 1/x prefactor in (f4.47|) 



which has its origin in the contribution of spin fluctuations. On the other hand, correlation 
functions that involve operators changing the total number of particles (e.g. the single 
particle Green function and pairing correlation functions) decay like exp[— est. (In x)~ 3 ^ 2 ], 
i.e. faster than any power law. This in particular means that the momentum distribution 
function rik and all its derivatives are continuous at kp, and there is only an essential 
singularity at kp. 

It is instructive to notice that the above result ( f4.47|) , obtained in the limit of weak 



Coulomb interactions, can equally be derived assuming strong long-range repulsion |42 
by expanding around the equilibrium configuration of an equally spaced crystal. We thus 
are lead to the rather remarkable conclusion that the long-distance behavior of correlation 
functions is independent of the strength of the Coulomb repulsion, provided the interaction is 
truly long-ranged. We note that an interpretation of Luttinger liquid properties in terms of 



a one-dimensional harmonic chain can also be given for short-range interactions p7],|28|,f4*6 



4- Persistent current 

The Luttinger model description can be used straightforwardly to obtain the current 
induced in a strictly one-dimensional ring threaded by a magnetic flux $. p7|,ff8| The ar- 



gument can in fact be made very simply: in the one-dimensional geometry, the vector field 
can be removed entirely from the Hamiltonian via a gauge transformation, which then leads 
to the boundary condition ijj(L) = exp(2TTi§/§o)ip(0) for the fermion field operator. Here 
L is the perimeter of the ring, and $o — hc/e. For spinless fermions, this is achieved by 
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replacing 



2$ 

H(x) - U(x) + — (4.48) 



in the bosonization formula, eq.( 4.19| ). The total J-dependent part of the Hamiltonian then 
becomes 

TT1I J\ 

Hj = —(J + 2$/%) 2 , (4.49) 

giving rise to a number current 

%dE uK ( 2$\ 

At equilibrium, J is chosen so as to minimize the energy. Given that at constant total particle 
number J can only change by two units, one easily sees that the equilibrium (persistent) 
current has periodicity <3>o> an d reaches is maximum value uK/L at $ = $0/2, giving rise 
to the familiar sawtooth curve for the current as a function of flux. 

For fermions with spin, as long as there is no spin gap (gi > 0), the above results can be 
taken over, with the simple replacement uK — > 2u p K p , the factor 2 coming from the spin 
degeneracy. Note in particular that the persistent current, an equilibrium property, is given 
by the same combination of parameters as the Drude weight in the conductivity. This is an 



illustration of Kohn's result |49| relating the Drude weight to the effect of a magnetic flux 
through a ring. 

In the case of negative g\, electrons can be transfered from the right to the left-going 
branch only by pairs, Consequently, the periodicity of the current and the ground state 
energy is doubled to 2$ , and the maximum current is equally doubled. This behavior has 
actually been found in numerical calculations. [f)0,[5l[ 



V. TRANSPORT IN A LUTTINGER LIQUID 



In the previous section we were concerned with equilibrium properties and correlation 
functions, in order to characterize the different phases possible in a one-dimensional system 
of interacting fermions. Here, we will investigate transport, in particular the dc conductiv- 
ity. Finite-frequency effects have also been investigated, and the reader is referred to the 
literature. [pTj^2| 



A. Conductance and conductivity 



To clarify some of the basic notions, let us first consider a Luttinger model in the presence 
of a weak space- and time-dependent external potential if. The interaction of the fermions 
with if is described by the extra term 

H ex t = ~e J dxp(x)(p(x,t) (5.1) 

in the total Hamiltonian. We will assume that the external field is slowly varying in space, 
so that in the particle-density operator p only products of either two right- or two left-going 
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fermion operators appear but no cross terms. Standard linear response theory tells us that 
the current induced by the potential is given by 



e 2 rt 



oo 



j(x, t) = -— I dt' I dx'D jp (x -x',t- t')if(x', t') , (5.2) 



where the retarded current-density correlation function is given by 

D jp (x,t) = -i9(t)([j(x,t),p(Q,Q)}) 
u K 

= — e-PQ(t) [S'( x - u p t) + 8\x + V)] • (5-3) 

7T 

The second line is the result for spin-1/2 electrons. For spinless fermions one has to make 
the replacement u p K p — > uK/2. 

Let us now first consider the situation where we adiabatically switch on a potential of 
frequency u and wavenumber q along the whole length of the system. From eq. (|5.3| ) one 
then straightforwardly obtains the q- and cu-dependent conductivity as 

/ \ 4e 2 Tyr i(uj + i5) . . 

<J\Sl^) = —u p K p - — — — — . 5.4 

k (<jj + ib) z — u A p q l 

In particular, the real part of the conductivity for constant applied field is 

2e 2 

<t(0,u) = — u p K p 5(u) , (5.5) 



in agreement with eq.( 4.37|) (where units with e 2 = h = 1 were used). 



Applying on the other hand a static field over a finite part of the sample, one obtains a 
current j = 2e 2 K p U / h, where U is the applied tension. The conductance thus is 

2e 2 

G = —K p , (5.6) 

and depends on K p only, not on u p . For the noninteracting case K p = 1 this is Landauer's 
well-known result. ]53| Note that interactions affect the value of the conductance. The 
conductance here is independent of the length over which the field is applied. Noting that 
in dimension d the conductance is related to the dc conductivity via G = L d ~ 2 a, a length- 
independent conductivity implies an infinite conductivity in one dimension, in agreement 
with eq.(|5.3|). The fact that u p does not appear in the expression for G can be understood 
noting that applying a static field over a finite (but large) part of the sample, one is essentially 
studying the wavenumber-dependent conductivity at strictly zero frequency, which from 
eq. (|5.15| ) is given by a(q, 0) = 2e 2 K p 5(q)/h, indeed independent of u p . On the other hand, 
applying a field of finite frequency over a finite length £, one can see that one measures the 
conductivity a(q — > 0, lj) if uj > u p jl. 

B. Edge states in the quantum Hall effect 

The Luttinger liquid picture has an interesting application to the physics of the fractional 
quantum Hall effect, as discovered and discussed by Wen. 0,^5j To see how this comes 
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FIG. 8. Quantum states in the lowest Landau level in the vicinity of the edge of a quantum hall 
device. The spatial variation of the confining potential is assumed to be slow on the scale of the 
magnetic length, so that the energies of the different quantum states are determined by the local 
value of the confining potential. For a circular device the angular momentum of a state increases 
proportionally to its distance r from the origin. 

about, consider the states available in the different Landau levels in the vicinity of the 
edge of the quantum Hall device, as shown in fig.[| It is clear that low-energy states 
only exist at the edge (the bulk quantum Hall state is well-known to be characterized by 
a finite excitation gap), and close to the Fermi energy (i.e. the edge) the states have a 
linear dispersion relation. This can be made particularly clear if one assumes a disk-shaped 
sample: the states have a wavefunction 

ss z k « (e ie ) k , 

with k increasing linearly with radial position. The angular momentum quantum number k 
thus plays a role very similar to linear momentum in the linear geometry we have assumed 
up to now. One thus can linearize the the dispersion in fig.|8| and obtains essentially the 
spinless model discussed in section [IV A|, the only difference being that here only right- 



going particles exist. This difference is the origin of the term chiral Luttinger liquid (in fact, 
the left going branch is to be found on the opposite edge of the device). Because there are no 
left-going particles (or at least they can be thought of as being at a macroscopic distance), 
there also is no right-left interaction, and consequently one expects the noninteracting value 
K = 1. Moreover all the left-going components of the fields have to be projected out, for 
example one has to replace — > <p + = (0 — 9)/2. 

However, straightforward adoption of this scheme leads to trouble: from the preceding 
subsection we know that K = 1 leads to a conductance (which in the present case is the 
Hall conductance) of G — e 2 /h, different from the well-known 

G = vj (5.7) 
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valid for a fractional quantum Hall state {y — 1/m is the filling factor). To repair this 
problem one makes the hypothesis that instead of eqs.( 4.18] ) and (|4.23|) one has 



( = _ v^^± ]{x) = u ^ U+ (x) , (5.8) 

7T ax 

where the subscripts indicate projection on right-going states. With these definitions fol- 
lowing the calculations of the previous subsection one now straightforwardly reproduces the 
correct result, eq.([0]). The appearance of the factors \/v in eq. (|5.8|) indicates that the ob- 
jects occupying the states in fig.g are not free electrons but rather strongly affected by the 
physics of the bulk of the samples. A more detailed derivation, starting from a Chern-Simons 
field theory for the bulk physics, has also been given by Wen. |55| 



Beyond reproducing the correct value of the Hall conductance, the above hypothesis 
leads to a number of interesting conclusions. Consider first the creation operator for a real 
electron (charge e) on the edge. Following the arguments of section |IV A] , because of eq. (|5.8| ), 



the bosonized version of the electron operator now must create a jump of of height vr/y^, 
rather than of height ir. This leads to 

ip+(x) « e - iMx)/ ^ , (5.9) 

x being the coordinate along the perimeter of the sample. Now, these operators obey the 
relation 

^+(x')^+(x) = e ±in/ ^ + (x)^ + (x') . (5.10) 

But the real electron is still a fermion, i.e. ip + must obey anticommutation relations. Thus 
m = 1/V has to be an odd integer. One thus reproduces one of the fundamental facts of the 
fractional quantum Hall effect. From eq. (|5.9| ) one also finds a decay of the single-electron 
Green function as 

G(x,t) oc — -r . (5.11) 

Another fundamental property of the quantum Hall state appears when one considers 
the fractionally charged elementary excitation of charge ev at the edge. A charge-e^ object 
is created by 

i/> +u (x) » e~ l ^ +{x) , (5.12) 

leading to a slow decay of the corresponding Green function, with exponent u, instead of 
1/u in eq.( |5.11| ). One now has the relation 

f/> +v (J)i/> +v (x) = e^j+AxW+vW ■ (5.13) 

i.e. the fractionally charged objects also obey fractional statistics. 

A single hypothesis, the insertion of the factors y/u in eq. (|5.8| ), thus reproduces two of 
the fundamental facts about the fractional quantum Hall effect! In addition one obtains 
results for the asymptotics of Green functions. 
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C. Single barrier 



The infinite conductivity in the ideally pure systems considered up to here is a natural but 
hardly realistic result: (almost) any realistic system will contain some form of inhomogeneity. 
This in general leads to a finite conductivity, and in one dimension one can anticipate even 
more dramatic effects: in a noninteracting system any form of disorder leads to localization 
of the single-particle eigenstates. [|57],[5!|] How this phenomenon occurs in interacting systems 
will be discussed in this and the following section. 

Following Kane and Fisher, [|5lJ consider first the case of a single inhomogeneity in an 
otherwise perfect one-dimensional system. The extra term in the Hamiltonian introduced 
by a localized potential v(x) is (for spinless fermions) 



Hbarrier = I dxiffi (x)l[> (x) . (5-14) 

Decomposing the product of fermion operators into right- and left-going parts, one has 

ifrty = ^ + tP + + ipl_ip_ + V>+0- + V'-'0+ • (5.15) 

In the bosonic representation, the first two terms are proportional to d x <p (see eq.( |4.18|) ), and 
therefore the corresponding contribution in eq.( |5.14| ) can in fact be eliminated by a simple 



unitary transformation of <p. These terms represent scattering with momentum transfer 
q <C 2kp, i.e. they do not transfer particles between kp and —hp and therefore do not affect 
the conductance in any noticeable way. On the other hand, the last two terms in eq.( |5.15|) 



represent scattering with \q\ m 2kp, i.e. from the + to the — branch and vice versa. These 
terms certainly are expected to affect the conductance, because they change the direction 
of propagation of the particles. The bosonic representation of these terms is 

Hbarrier = Y——?- COS 20(0) , (5.16) 
TTOi 

where the potential V(x) is assumed to be centered at x — 0. For this reason, only the value 
of the at x = intervenes in eq.( 5.16j ). 



One now can integrate out all the degrees of freedom away from x = 0, to obtain an 
effective action implying only the time-dependence of 0(0). Then a renormalization group 
equation for V = V(2kp) can be found as 

^ = (1-K)V, (5.17) 

where E = E Q e~ e , E is the original cutoff, and E is the renormalized cutoff. 
From eq . ( 5 . 1 7|) it is clear that there are three regimes: 



1. For K > 1 one has V(£ — > oo) = 0, i.e. as far as the low-energy physics is concerned, 
the system behaves like one without the barrier. In particular, the low-temperature 
conductance takes the "pure" value G = e 2 K/h, with corrections of order T 2<yK ~ l \ 
We note that in this case superconducting fluctuations dominate, and the prefect 
transmission through the barrier can be taken as a manifestation of superconductivity 
in the one-dimensional system. 
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2. For the noninteracting case K = 1, V is invariant, and one thus has partial transmission 
and a non-universal conductance depending on V. 



3. For if < 1 7(1) scales to infinity. Though the perturbative calculation does not 
provide any direct way to treat this case, it is physically clear that the transmission 
and therefore the conductance should vanish. 

Note that the non-interacting case is marginal, separating the regions of perfect and zero 
transmission. These results are very similar to earlier ones by Luther and Peschel |BIJ who 
treat disorder in lowest-order perturbation theory. 

The case of K < 1 can be further analyzed considering the case of two finite Luttinger 
liquids coupled by a weak tunneling barrier, as would be appropriate for a strong local 
potential. The barrier Hamiltonian then is 

H bamer = #1(0)^(0) + Vl(0)Vi(0)] « — cos2#(0) . (5.18) 

Here 2 are the field operators to the left and to the right of the barrier. The operators 
have to satisfy the fixed boundary condition [x — 0) = 0, different from the periodic 
boundary conditions we have used so far. Noting that the ip i can be decomposed into left— 
and right-going parts as ip i = ip^ + ip^, and using eq.( |4.19| this can be achieved by imposing 



the fixed boundary condition = 0) = n/2 on the boson field. f6l"|j62| 

One can now proceed in complete analogy to the weak-V case to obtain the renormal- 
ization group equation 

^ = (1-1/K)t . (5.19) 

Again, there are three different regimes: (i) for K > 1 now t(£ — > oo) — > oo, i.e. the tunneling 
amplitude becomes very big. This can be interpreted as indicating perfect transmission, e.g. 
G = e 2 K/h; (ii) the case K = 1 remains marginal, leading to a t-dependent conductance; 
(iii) for K < 1 t scales to zero, there thus is no transmission, and G — 0. 

The results obtained in the two limiting cases of small V (weak scattering) and of small t 
(weak tunneling) are clearly compatible: e.g. for K < 1 V becomes large, i.e. at sufficiently 
low energies one expects essentially a tunneling type behavior, and then from eq.( |5.19| ) the 
tunneling amplitude actually does scale to zero, giving zero conductance in the low-energy 
(or low-temperature) limit. For K > 1 a similar compatibility of the two limiting cases is 
found. The global behavior can be represented by the "phase diagram" in fig.[|. 

For electrons with spin but spin-independent interactions, results are very similar: the 
separation between zero and perfect transmission is at K p = 1, with K p = 1 again the 
marginal case. In the transmitting region the conductance is G = 2K p e 2 /h. 

These considerations can be generalized to the case of two barriers. p9]J5B[ | In particular, 
assuming that there are two identical, weakly scattering barrier at ±d, the effective scattering 
potential becomes V e ff(q) = 2v(q) cos(qd/2). Though in general this is non-zero when V(q) 
is non-zero, for particular values of kp, so that cos(&pd) = this potential vanishes, giving 
rise to perfect transmission even for K < 1. This resonant scattering condition corresponds 
to an average particle number between the two barriers of the form v + 1/2, with integer v, 
i.e. the "island" between the two barriers is in a degenerate state. If interactions between 
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FIG. 9. "Phase diagram" of a localized inhomogeneity in a spinless Luttinger liquid, character- 
ized by an exponent parameter K, according to [p9[] . The scaling trajectories calculated for weak 
V or t are indicated by arrows. It is clearly plausible to assume direct scaling from weak to strong 
coupling in the whole range of K. 



the electrons in the island are included, one can recover the physics of the Coulomb blockade. 

urn 



For the chiral Luttinger liquid discussed in the preceding section backscattering events a 
priori seem to be excluded because all the particles are moving in the same direction. In that 
sense the chiral Luttinger liquid can be considered as "perfect". However, if the quantum 
hall device has a constriction that brings the two edges close to each other, scattering from 
one edge to the other becomes possible and is the equivalent of backscattering. Then similar 
considerations as made for the single-impurity case are possible ||64j| , and in particular 
the crossover function describing the conductance through a resonance as a function of 
temperature has been obtained. pS| 



D. Random potentials and localization 

The discussion of the previous section was concerned with the effect of a single impurity, 
weak or strong. Clearly, in that case the effects of coherent scattering from many impurities 
which typically give rise to Anderson localization are absent. We now turn to this more 



complicated case which had been studied in fact well before the single impurity work. [|6r^-p9 
In the absence of electron-electron interactions, localization effects can be discussed in 
the framework of a scaling theory. ]53 Under the assumption that at some short length 



scale one has elastic scattering of electrons off impurities, this theory leads to the following 
/3-function for the variation of the conductance with linear dimension L: 

where a is a constant and d the spatial dimensionality. In particular in one dimension this 
leads to a conductance decaying exponentially with the length of the system, exhibiting 
clearly the localized character of the all single-electron states (a fact first shown by Mott 
]57| and studied in great detail since |]7D|,|7TJ). 



In an interacting one-dimensional system (as described by the Luttinger liquid picture 
of the previous section) now a number of questions arise: what is the influence of disorder on 
the phase diagram obtained previously? What are the transport properties? Can one have 
true superconductivity in one dimension, i.e. infinite conductivity in a disordered system? 
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To answer these question we discuss below the generalizations necessary to include disorder 
in our previous picture. 

We start by the standard term in the Hamiltonian describing the coupling of a random 
potential to the electron density 



Himp 



J2 I dx V(x - Ri)p(x) , (5.21) 
Ri J 



where the Ri are the random positions of impurity atoms, each acting with a potential V 
on the electrons. In one dimension one can distinguish two types of processes: (i) forward 
scattering, where the scattered particle remains in the vicinity of its Fermi point. As in the 
single-impurity case, this leads to a term proportional to <90 p and can be absorbed by a 
simple redefinition of the 4> p field. The physical effects are minor, and in particular the dc 
conductivity remains infinite, (ii) backward scattering where an electron is scattered from 
hp to — kp or vice versa. For small impurity density this can be represented by a complex 
field £ with Gaussian distribution of width = UiV{q = 2kp) 2 : 

H b = Y, [dx [e(x)4 CT (x)^ LCT (x) + h.c] (5.22) 

This term has has dramatic effects and in particular leads to Anderson localization in the 
noninteracting case. |72 



From a perturbative expansion in the disorder one now obtains a set of coupled renor- 



malization group equations |69j : 



dl 2u a p 

dK a 1 tin 2 \ 



dy 

di 
dV 

~dJ 



--{V + y*)K$ 
2{\-K a )y-V 

— K p — K a — y)V (5.23) 



where V = 2D^a/ '{^u 2 a ) {u a /up) Kp is the dimensionless disorder, y = gi±/(nu a ) is the di- 
mensionless backscattering amplitude, and the K v are defined in eq. (|4.32| ). These equations 



are valid for arbitrary K v (the usual strength of bosonization), but to lowest order in T> and 

y- 

As a first application of eqs. (|5.23|) one can determine the effect of the random potential 
on the "phase diagram", as represented in fig.[7]. In fact, there are three different regimes: 

1. for K p > 2 and g±± sufficiently positive the fixed point is V*, y* =0, K* > 2. Because 
the effective random potential vanishes this is a delocalized region, characterized as 
in the pure case by the absence of a gap in the spin excitations and dominant TS 
fluctuations. 

2. For Kp > 3 and g\± small or negative one has T>* = 0, y — > — oo, K* > 3. Again, 
this is a delocalized region, but now because y — > — oo there is a spin gap and one has 
predominant SS fluctuations. 
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3. In all other cases one has T> — > oo,y — > — oo. This corresponds to a localized regime. 
For small if p the bosonized Hamiltonian in this regime is that of a charge density 
wave in a weak random potential with small quantum fluctuations parameterized by 
K p . This region can therefore be identified as a weakly pinned CDW, also called a 
"charge density glass" (CDG). The transition from the CDG to the SS region then 
can be seen as depinning of the CDW by quantum fluctuations. 

One should notice that the CDG is a nonmagnetic spin singlet, representing approximately a 
situation where localized single-particle states are doubly occupied. Though this is accept- 
able for attractive or possibly weakly repulsive interactions, for strong short-range repulsion 
single occupancy of localized states seems to be more likely. One then has a spin in each 
localized state, giving probably rise to a localized antiferromagnet with random exchange 
(RAF). A detailed theory of the relative stability of the two states is currently missing and 
would certainly at least require higher-order perturbative treatment. The boundaries of the 
different regimes can be determined in many cases from eqs.( |5.23| ), and the resulting phase 



diagram is shown in figJTO 
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K„ 



SS 



FIG. 10. Phase diagram of a Luttinger liquid in the presence of a weak random potential 
(T> = 0.05). The full lines represent results obtained directly from the scaling equations ( 5.23| ), the 
dashed lines are qualitative interpolations. The dotted lines are the phase boundaries in the limit 
V -» 0. 

The localization length for small disorder can be obtained from standard scaling argu- 
ments: suppose that a system with some fixed disorder T>q has a localization length £o- Then 
in the general case one has £(T>) = £q^ v ° , ' d \ where £(Vo,T)) is the time it takes for the 
"bare" disorder D to scale up to T>q. From this reasoning one finds, for the case without a 
spin gap (pi > 0) and weak disorder 

£{V) oc (l/X?) 1 ^ 2 -^) . (5.24) 

Note that for K p > 1, i.e. superconducting fluctuations predominating in the pure case £ is 
greater than the mean free path A oc 1/T>, i.e. there is a kind of diffusive regime, contrary 
to the noninteracting case. On the other hand, for K p < 1 one has £ < A. In the vicinity of 
the TS-RAF boundary one has 



£{V) oc exp 



V-y(K p -2) i 



(5.25) 



39 



The analogous results for the case with a spin gap (gi < 0) are 



i{V) oc (1/D)V(3-^) (5.26) 
iiV) oc exp { , 2n I (5.27) 

There are two points to be noted about this result: (i) for K p = one has £ oc "D -1 / 3 , 



which is the same result as that found for the pinning length of a classical CDW. 173 



n 



the results (|5.24|) , ( |5.25| ) and ( |5.26|) , eq. (|5.27|) are qualitatively different, both in the vicinity 



of the phase boundaries and in the localized states. The transitions are thus in different 
universality classes, and this strongly supports the idea that the localized phases reached 
through the transition are themselves different (RAF or CDG). 

The temperature dependence of the dc conductivity can be obtained noting that at finite 
temperature there are no coherent effects on length scales larger than vf/T. One therefore 
stops renormalization at e r = Vf/{olT). As long as D remains weak one can still use the 
Born approximation to obtain 

= a ^ , (5.28) 

where gq = e 2 v F /2irhD^ is the lowest order conductivity. In the delocalized phases one 
then finds a conductivity diverging as cr(T) ~ T~ 1-7 where 7 = K* — 2 in the TS case 
and 7 = K* — 3 in the SS case. On the phase boundaries one has universally o ~ 1/T. 
In the localized region V diverges at low temperatures, and a perturbative calculation thus 
becomes meaningless. However, the conducting-localized crossover can still be studied at 
not too low temperatures. |){| In particular, the high-temperature conductivity is found 



to vary as a ~ T 1 Kp . This is the perturbative result first found by Luther and Peschel 
IBTjfl and also reproduced by the single-impurity calculations. J59J The high-temperature 



behavior thus can be understood in terms of scattering off the individual impurities. On the 
other hand, at lower temperatures one necessarily comes into the region where T> increases 
sharply. This has its origin in coherent scattering from many impurities and ultimately gives 
rise to localization. 

One can finally notice the effects of different types of interactions on localization. Roughly 
speaking, for forward scattering repulsion (g 2 > 0) enhances localization whereas attraction 
weakens it. In particular, strong attraction leads to vanishing effective random potential. 
The delocalized state then can be considered to be a true superconductor in the sense that 
there is infinite conductivity even in an impure system. The effect of backward scattering 
interactions is opposite to that of forward interactions. 



VI. THE BETHE ANSATZ: A PEDESTRIAN INTRODUCTION 

The methods and results described in the preceding two sections were to a large extent 
based on a continuum description of the physics of one-dimensional electrons, and often 
relied on perturbation theory. Complementary insight can be gained by considering models 
defined on a lattice, and in particular some exactly solvable models. The physically most 
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important cases are the Heisenberg spin chain (the subject of Bethe's original work) and 
its generalization for itinerant fermions, the Hubbard model. In the following, I will dis- 
cuss the exact solution of the Heisenberg model in some detail, in the hope to clarify the 
essential features. Subsequently, the necessary generalizations for the Hubbard model will 
be indicated. 

The exact solutions give exact energies of the ground state and all the excited states 
in terms of the solution of a system of coupled nonlinear equations. On the other hand, 
the corresponding wavefunctions have a form so complicated that the explicit calculation of 
matrix elements, correlation functions and other physical quantities has remained impossible 
so far. In the subsequent sections I shall describe how the knowledge of the energy spectrum 
obtained from exact solutions can be combined with the results of the preceding two chapters 
to obtain a rather detailed picture of the low-energy properties, in particular of correlation 
functions. 

It should be noticed that Bethe's method has been successfully applied to a number of 



other interesting cases like the Kondo model [74,75] and also is by no means restricted to 



lattice models. [75-7§| For more general and detailed discussions, the reader is referred to 
the literature. 



U-0 



A. The Heisenberg spin chain 



1. Bethe's solution 



We will write the Hamiltonian of the Heisenberg spin chain of L sites in the form 



H — E(^ X ^i+i + SfSf +1 + AS'fS', 



z N 

1+1/ 



i=l 
L r 



i=l lZ 



(6.1) 



Here Si = (Sf, Sf , Sf) is a spin-1/2 operator acting on site i, S is an anisotropy parameter 
that allows one to treat the antiferromagnetic (A = 1), the ferromagnetic (A = —1), and 
general anisotropic casesQ, and periodic boundary conditions imply Sl+i = Si- We notice 
that the Hamiltonian conserves the z- component of the total spin (for |A| = 1 total spin is 
also conserved), and consequently one can write any eigenfunction of a state with N down 
and L — N up spins as 



E f(x 1 ,x 2 ,...,x N )\x 1 ,x 2 ,...,x N ) 



(6.2) 



Xl<X2<--.<X N 



2 In the Hamiltonian, eq.( |6.1| ), an overall sign opposite to that of some of the original literature 
Hp4| is used 
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where the Xi are the positions of the down spins. Applying the Hamiltonian to this state, 
one obtains a "Schrodinger equation" for the coefficients /. In particular, if none of the x^ 
are adjacent, one has 

1 N 

*/({**}) = 5 £*</({*<}) > (6-3) 
A i=i 

where e = E + A(N — L/4), E is the eigenenergy, and the finite-difference operator 8^ is 
defined by 

8if(. ....,;....) /( r,- 1....) • /( r, • 1....) . (6.4) 

If two down spins are adjacent, e.g. x k+ i = x k + 1, displacement of the spin is partly 
suppressed, and one finds instead of eq.(|6.3|) 

1 N 

sf(..., x k ,x k + l,...) = - 8if(. . . ,x k ,x k + l,. . .) 

+ 2 [/(■ • • > x fe _ !> x fc + !> • • •) + /(• • • > x k, %k + 2, ■ ■ •)] 

+ Af(...,x k ,x k + 1,...) (6.5) 

It should be quite clear how this generalizes to the case of more than two Xj's adjacent: 
whenever an down spin has an down neighbor, the term representing displacement to that 
site disappears, and one finds an "interaction term" (A) instead. 

Bethe's solution starts by the observation that any linear combination of plane waves of 
the form 



f(x u ...,x N ) = ^Apexp 
p 



N 

i £ kpjXj 



(6.6) 



is a solution of the "free" eq.( |6.3|) for arbitrary coefficients Ap. Here the summation is over 
all permutations P of the iV different wavenumbers kj, i.e. in eq.( |6.6D the momenta are 
permuted amongst all the particles (down spins). 

The essential point in Bethe's solution is now to determine the coefficients Ap so that 
eqs.(|6.3D and (p.5|) become identical, i.e. one requires 



2A/(. . .,x k ,x k + l, . . .) = 

/(. .. ,x k + l,x k + 1, ...) + /(... ,x k ,x k , .. .) . (6.7) 

Inserting this into eq. (|6.5| ) one sees that one recovers eq.( |6.3|) which is automatically solved 
by the plane wave form (|6.6|) . One should notice in eq. (|6.7|) that coefficients with two 
equal arguments like /(. . . ,x k ,x k , . . .) have no direct physical meaning but are of course 
well-defined mathematically from eq. ( |6.6D . 

To see how this works let us first consider the two-particle case. Then there are only 
two permutations to sum over in eq. (|6.6|) , and eq. fl6.7|) becomes (after canceling common 
factors) 

2A(A e lk2 + A ie lkl ) = (Aq + A 1 )(l + e l(kl+k2) ) . (6.8) 
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This then gives 

A\ 1 + e*( fcl+fc2 ) — 2Ae ifc2 

T = - i + e i(*i +fc3 )-2Ae^ = " exp H©(^, h)} , (6.9) 

where the second equation defines the phase shift G. Note that this means that the plane- 
wave form of the solution, eq.( |6.6|) , remains valid even when two down spins are at nearest- 
neighbor sites where they interact. For longer range interactions this would not in general 
be the case. 

We now consider eq.( |6.7| ) for general N. It is then convenient to associate with each 
permutation P another permutation P' which differs from P only by the exchange of two 
adjacent elements: P'k = Pk + 1, P'k + 1 = Pk. One then has 



/(..., x k + m, x k + n, ...) = Yl' ex P 



Q i{k Pk +k p i k )x k 



I k P3 X 3 

x (Ape^ mkpk+nk P'" ) + A P ,e i{mkp ' k+nkpk) ^ , (6.10) 

where in eq.( |6.7| ) the cases m, n = 0, 1 are relevant, and the summation is over half of the 
permutations (for example those with Pk < P'k), the other permutations being included 
explicitly by the term proportional to Api. Given that eq.( |6.7|) is supposed to be valid for 
any set {xi}, the coefficients of the plane wave factors in the left and right hand sides have 
to be equal, e.g. 

2A(A P e ik P'>° + A P ,e ikpk ) = {A P + A P ,){1 + e i{kpk+kp ' k) ) . (6.11) 
This is an obvious generalization of eq. (|6.9|) and leads to 

A p i -f e i( k Pk+k P / k ) _ 2Ae ikpk 

AT = " l + e ^+W-2Ae^ = "° pHe( * w,W) ' ^ 

This now fixes all the coefficients in eq.( |6.6| ) unambiguously, up to an overall normalization, 
because every permutation can be build up from a sequence of elementary permutation of two 
adjacent elements. To each such sequence corresponds a product of factors like eq. (|6.12|) , one 



for each elementary permutation, and moreover this product is independent of the sequence 
of individual permutations, e.g. 

(123) -> (213) -> (231) -»■ (321) 
and (123) -> (132) -> (312) -> (321) (6.13) 

lead to the same factor. One should also notice that eq.( |6.12j ) was derived without any 
assumption on the other x^s, and consequently situations with three or more adjacent down 
spins are equally covered and the generalizations of eq. (|6.5|) to these cases are therefore also 
fulfilled. 

It remains to determine the allowed /c-values. This follows from the periodic boundary 
condition 

f(l,x 2 , • • • ,x N ) = f(x 2 , ■ ■ ■ ,x N ,L + 1) . (6.14) 
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Inserting into eq.( |6.6|) and noting that eq.( |6.14|) has to be satisfied for all (X2, . . . ,%n), one 
obtains the condition 



Ap/A P ,e ik ^ L = 1 



(6.15) 



which must be satisfied for every permutation P. Here P' is a permutation obtained from 
P by a "right shift" of all elements: 



P' = (PN, PI, P2, . . . , PN - 1) . 



(6.16) 



The ratio of coefficients in eq.( |6.15| ) now can be calculated by repeatedly permuting PN in 
eq. (|6.16|) to the right and using eq. fl6.12|) . The resulting equation is 



N 



-l) N - 1 exp\iY l &(k j ,k l )\e iklL = l , 

3=1 



(6.17) 



with defined in eq. fl6.12p . 



2. The Heisenberg antiferromagnet 

We now turn to the particularly interesting and physically relevant case of the Heisenberg 
antiferromagnet (A = 1). In this case, the fundamental eq. (|6.17|) can be transformed into a 
somewhat simpler form making the transformation of variables 

A = -icot(fc/2)^e* fc = ^ . (6.18) 

Equation (|6.17 ) then becomes 



L at 



'2A,-A ^Xi-Xj 



\2Xi + iJ jJ{Xi-Xj + i 

and similarly one finds the total momentum 

N N 

P = J2kj = J2( 2 arctan(2A j ) - tt) (6.20) 
i=j j=i 

and the energy 

N N 2 

l2 



E = L/A + £(cos % _ 1) = L/4 - X: T-7T2 

3=1 3=1 1 + 4A J 



(6.21) 



The z-component of the total spin of the states obtained from eq.( |6.19| ) is obviously L/2 — N. 
In fact it can be shown that this is also the total spin of these states, and that these states, 
together with those obtained by repeated application of the lowering operator of the total 
spin, form the complete set of eigenstates. 

If all the A's are real, eq. (|6.19|) can be simplified by taking the logarithm. One then 
obtains 

N 

2L arctan(2Ai) = 2vr J x + 2 ^ arctan(A ; - X A ) , (6.22) 

3=1 
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where the J\ are integers if L — N is odd and half-odd integers (of the form 1/2, 3/2, ...) is 
L — N is even and are restricted to 



| J, | < (L-N + l)/2 . (6.23) 

In the ground state one has N = L/2, and then all the Ji are uniquely determined. In the 
thermodynamic limit, the allowed A's are very close to each other, and their distribution 
then can be determined by a linear integral equation. The ground state energy then is 
E = L(l/4-ln2). 

Excited states can be obtained by taking out one of the down spins. One then has a 
state of total spin 1. One now has to choose N = L/2 — 1, subject to the constraint ( |6.23D , 



i.e. one has to take two (half-) integers out of the sequence —L/A...L/A. This is thus a 
two-parameter family of states, which, in the thermodynamic limit has energy 

7T 

E(Pi,P2) = e(pi) +e{P2) , with e(p) = -sinp (6.24) 
and momentum P = p\ + p 2 . At fixed total momentum P there is thus a continuum 



of allowed states, as shown in figJTl]. This energy-momentum relation suggests that the 



excitation is actually the combination of two elementary excitations which don't interact, 




FIG. 11. The continuum of low-lying triplet and singlet states. The lower and upper limit of 
the continuum are (tt/2) | sin k\ and vr | sin(£;/2)|, respectively. 

so that there energies and momenta just add up. This interpretation finds further support 
in the fact that excitations with total spin exist that have the same energetics as the spin 
1 excitations (these correspond to the appearance of one pair of complex conjugate A's in 
eq. ( |6.19| )). These singlet and triplet excitations thus can in fact be seen as excitations of 
pairs of spin 1/2 objects, called "spinons", which do not interact but form states of total 
spin or 1. Spinons can be visualized as kinks in an antiferromagnetic order parameter, as 



shown in fig.|12| In an isotropic Heisenberg chain, these individual spinons get delocalized 
into plane wave states due to the exchange interaction. For an even number of sites the total 
spin is always integer, so that the spinons always have to be excited pairwise. However, once 
created, they behave essentially as noninteracting objects. 

In the Heisenberg model, the spinons are only noninteracting in the thermodynamic limit 
L — > oo. However, in a modified model with 1/r 2 interaction, this separation is true for all 
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(a) 



(b) 



FIG. 12. Two-spinon configurations in an antiferromagnetic chain with m z = (a) and m z = 1 



(b). 



states even in a finite system [^,^]. One can further show that individual spinons behave as 
semions, i.e. their statistical properties are in many respects intermediate between fermions 



and bosons IB? . 



3. The Jordan- Wigner transformation and spinless fermions 



The spin model discussed above can be transformed into a model of spinless fermions, 
noting that and S~ anticommute. The Jordan- Wigner transformation then relates 
spin to fermion operators (a i; aj) via 



a\ exp 



i-1 



1 71 



a j a j 



QZ 



(6.25) 



Presence or absence of a fermion now represent an up or down spin, and the exponential 
factor insures that spin operators on different sites commute, whereas fermionic operators of 
course anticommute. The transformation can now be used to rewrite the spin Hamiltonian 



(ITT) in terms of fermions as 



ri l l 



(6.26) 



The spin model thus is transformed into a fermionic model, with the "spin-flip" terms giving 
rise to motion of the fermions, whereas the S z -S z interaction gives rise to a fermion-fermion 
interaction between adjacent sites. 

The fermionic analogy now allows one to study correlation functions of the spin chain 
using the bosonization method developed in section [IV A| . In particular, the spin-spin cor- 
relations for the isotropic antiferromagnet decays as |BS| 



(^-^>«(-l)l^l 



In 1 / 2 


i - 






i - 


-3 





(6.27) 



where the logarithmic corrections term comes from contributions of a marginally irrelevant 
umklapp interaction term p8| , |37| , |39| . The alternation in eq. ( |6.27| ) indicates the expected 
tendency towards antiferromagnetic order, but correlations do decay (with a rather slow 
power law), and there thus is no long-range order, as to be expected in one dimension. 
The analysis can be extended to give the full dynamic structure factor p0 |, and recent 
experiments on KCuF 3 confirm these results. |H[ 
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For A < 1, a similar law holds, but with a A-dependent exponent and no logarithmic 
term. On the other hand, for A > 1 the spins are preferentially aligned along the ^-direction, 
and one then has a long-range ordered ground state of the Ising type. There thus is a phase 
transition exactly at the isotropic point A = 1. In the fermionic language, this corresponds 



to a metal-insulator transition 92 



B. The Hubbard model 

The Hubbard model is the prototypical model used for the description of correlated 
fermions in a large variety of circumstances, ranging from high-T c superconductors to heavy 
fermion compounds and organic conductors. In spite of its apparent simplicity, there is still 
no general solution, or even a consensus on its fundamental properties. Notable exceptions 
are the cases of one and infinite dimensions P^J54]] . In particular, in one dimension an exact 
solution is available. [p5 

The Hamiltonian in one dimension has the well-known form 

H = -tJ2( a ls a i+l,s + 4+l,s a i,s) + U J2 n iA n i,l > ( 6 - 28 ) 
i,s i 

where rij )S is the number operator for fermions of spin s, and the sum is over the L sites of 
a one-dimensional chain with periodic boundary conditions. 



The model has two global SU(2) symmetries p6|-|98|| : the first is the well-known spin 
rotation invariance and the second is particular to the Hubbard model and relates sectors of 
different particle numbers. The total symmetry thus is SU(2) x SU (2) ~ SO (4). One should 
notice that more complicated interactions, e.g. involving further neighbors, will conserve 
the spin rotation invariance but in general not the "charge" SU(2) invariance. Rather, this 
second symmetry will become the standard global £7(1) invariance associated with particle 
number conservation. 



1. The exact solution 

The exact solution of the one-dimensional Hubbard model has been found by Lieb and 



Wu |[95| . To obtain the wavefunctions of the Hubbard model, one writes a general iV-particle 
state as 

N 

\ F ) = E • • • E E • • • E r - v ) 1 1 < jSi \vac) . (6.29) 

XI X N Si s N i = l 

The simplest case is the two-particle problem. Then in the two parts of configuration space 
X2 > x\ (region I) and x\ > X2 (region II) the wavefunction is just a product of two plane 
waves, and the only nontrivial effects occur for x\ = x^- The full wavefunction F then is 

F SliS2 = e^ +k ^[9(x 2 - *0C 2 + - *2)CJ > ( 6 - 3 °) 
where the factors £ depend only on the spin quantum numbers of the two particles, and 9 is 
the usual step function, with 9(0) = 1/2. In order to satisfy the Schrodinger equation, these 
coefficients have to obey the equation 

£si,S 2 = E ^SlS2,ttt2^tlt 2 • (6.31) 

ti,h 
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The S-matrix has the form 

_ t(sin ki - sin k 2 ) + iU P a p a5 , , 

t{sm k\ — sm K2J + iU 

where Pa/3,75 is the operator permuting two spins. This operator acts on the spin part of the 
two-particle Hilbert space. 

For N particles, similar to the Heisenberg case, one has matching conditions whenever 
two particles cross. The major complication in the Hubbard model comes from the fact that 
these matching conditions obviously involve the spin of the two particles and therefore involve 
the S-matrix, eq. QB.32p , rather than just complex phase factors (compare with eq.(|6.9|)). In 
order to fulfill compatibility conditions for different paths in configuration space like those 
of eq. (|6.13|) one has has to have 



S 23 S 13 S 12 = S 12 S 13 S 23 . (6.33) 



These are the famous Yang-Baxter equations [[TTj] which have to be satisfied for a system 
to be solvable by Bethe Ansatz. One can verify that these equations are indeed satisfied 
by the S-matrix of the Hubbard model, eq. (|6.32|) . The subsequent analysis is based on a 



"generalized Bethe ansatz" [fTSLfTT]] . A detailed description of the methods used can be found 



in the specialized literature [80,81], and a very detailed derivation is given in Gaudin's book 



Imposing periodic boundary conditions, the allowed values of kj are obtained from the 
solution of the coupled set of nonlinear equations 



e ik 3 L = TT e I " a > ] (6.34) 



3=1 




flei ^-J™^ ) = -Ue[ ) , (6.35) 



Here N is the total number of electrons, M is the number of down-spin electrons (M < N/2), 
and e(x) = (x + i)/(x — i). The X a are parameters characterizing the spin dynamics. We 
note that in general, both the k/s and the A's are allowed to be complex. The energy and 
momentum of a state are 

N N 

E = -2tJ2 cos kj , P = Y. k J ■ ( 6 - 36 ) 



The determination of all the solutions of eqs. ( |6.34| , |6.35|) is not easy. It has recently been 



shown that under certain assumptions these equations do indeed give all the "lowest weight" 
(with respect to SU(2) x SU(2)) eigenstates of the Hubbard model. The complete set of 
eigenstates then is obtained acting repeatedly with raising operators on the Bethe ansatz 
states. [|9J 

In the ground state all the fc's and A's are real. Numerical results for the ground state 
energy as a function of particle density and U have been given by Shiba [ |10U|| . The excited 
states are discussed in some detail in the literature [|100| - |103| , ^6| . I here just briefly summarize 
the different types: 
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1. "4/cf" singlet states. The excitation energy goes to zero at q = and q = Akp. These 
excitations form a two-parameter continuum and only exist away from half-filling. 



2. u 2kp" singlet and triplet states. These are degenerate two-parameter families of solu- 
tions, similar to the Heisenberg chain. The energy goes to zero for q = and q = 2kp. 

3. There are also states with complex k. The energy of these states is proportional to U 
for large U : t, i.e. some sites are doubly occupied ("upper Hubbard band"). 

In addition, there are of course states with added particles or holes. 

For U — *■ oo the ratio sin kj/U in eqs.( |6.34] , |6.35[ ) clearly vanishes, however there is no 
restriction on the A's. Introducing the scaled variables A a = 2X a /U, the "spin equation" 
( |6.35| ) becomes 

2A Q + t \ N _ ^ A a -A f3 + t 

2A a —i) ^Av-Ap-i ' 1 ' ' 

which are the well-known Bethe ansatz equations for the spin-1/2 Heisenberg chain, 
eq. (|6.19|) , i.e. the spin wavefunction of iV particles is just that of an iV-site (not L-site!) 



Heisenberg chain, even when there is less than one particle per site. 



2. Luttinger liquid parameters 

The above results provide a complete picture of the energy spectrum of the one- 
dimensional Hubbard model (and also of its thermodynamics ||104|| ). On the other hand, 
the wavefunctions obtained are prohibitively complicated, and no practical scheme for the 
evaluation of expectation values or correlation functions has been found. Of course, in a 
weakly interacting system the coefficients K p and u v can be determined perturbatively. For 
example, for the Hubbard model one finds 

K p = l- U/(nv F ) + ... , (6.38) 

where Vp = 2t sin(7ra/2) is the Fermi velocity for n particles per site. For larger U higher 
operators appear in the continuum Hamiltonian ( 4.3U| ), e.g. higher derivatives of the fields 



or cosines of multiples of \/80 CT . These operators are irrelevant, i.e. they renormalize to zero 
and do not qualitatively change the long-distance properties, but they do lead to nontrivial 
corrections to the coefficients u u , K p . In principle these corrections can be treated order by 
order in perturbation theory. However, this approach is obviously unpractical for large U, 
and moreover it is quite possible that perturbation theory is not convergent. To obtain the 
physical properties for arbitrary U a different approach [|105|1 can be used which will now be 
explained. 

I note two points: (i) in the small-?/ perturbative regime, interactions renormalize to 



the weak-coupling fixed point g\ = 0, K* = 1; (ii) the exact solution |)5| does not show 
any singular behavior at nonzero U, i.e. large U and small U are the same phase of the 
model, so that the long-range behavior even of the large U case is determined by the fixed 
point g\ = 0. Thus, the low energy properties of the model are still determined by the three 
parameters u PjCT and K p . 



49 



0.0 0.2 0.4 0.6 0.8 1.0 
FIG. 13. The charge and spin velocities u p (full line) and u a (dashed line) for the Hubbard 
model, as a function of the band filling for different values of U/t: for u a U/t = 1, 2, 4, 8, 16 from 
top to bottom, for u p U/t = 16, 8, 4, 2, 1 from top to bottom in the left part of the figure. 

The velocities u P)Cr can be obtained from the long wavelength limit of the excitations 
discussed above. In particular, the Akp component of the charge correlation function only 
involves the (f) p field, and it is therefore natural to obtain u p from the q — > limit of the 
"4/cp" singlets. Similarly, u a is obtained from the "2fcp" excitations. In the thermodynamic 
limit the corresponding excitation energies are easily found from the numerical solution of 
a linear integral equation. [TUT] Results are shown in fig.[T3| for various values of U/t. Note 



that for U = one has u p = u a = 2t sin(7m/2), whereas for U — > oo u p = 2tsin(7rn) 
= (2irt 2 /U)(l — sin(27rn)/(27rn)). In the noninteracting case u a oc n for small n, but for 

any positive U u a oc n 2 . 

To obtain the parameter K p from the exact solution note that the gradient of the phase 

field <p p is proportional to the particle density, and in particular a constant slope of <f) p 

represents a change of total particle number. Consequently, the coefficient u p /K p in eq. 



(4.31) is proportional to the variation of the ground state energy Eq with particle number: 



1 d 2 E (n) vr u p 1 

L On 2 2K p n 2 K ' { ' 

Equation ( |6.39| ) now allows the direct determination of K p : E (n) can be obtained solving 
(numerically) Lieb and Wu's |95j integral equation, and u p is already known. The results for 



K p as a function of particle density are shown in fig.[L4] for different values of U / 1. For small 
U one finds in all greement with the perturbative expression, eq. (|6.38| ), whereas 

for large U K p — > 1/2. The limiting behavior for large U can be understood noting that 
for U = oo the charge dynamics of the system can be described by noninteracting spinless 
fermions (the hard-core constraint then is satisfied by the Pauli principle) with hp replaced 
by 2kp. Consequently one finds a contribution proportional to cos(4kpx)x~ 2 in the density- 
density correlation function, which from eq. (|4.42|) implies K p = 1/2. One then finds an 
asymptotic decay like cos(2/ci?x)a;~ 3//2 ln 1 ^ 2 (x) for the spin-spin correlations, eq. (|4.43 ), and 



an exponent a = 1/8 in the momentum distribution function. The result a = 1/8 has also 



been found by Anderson and Ren, ||106|| and by Parola and Sorella | 107|| . Ogata and Shiba's 
numerical results | 108 | are quite close to these exact values. 

We note that in the whole parameter region, as long as the interaction is repulsive 
one always has K p < 1, which means that magnetic fluctuations are enhanced over the 
noninteracting case. On the other hand, superconducting pairing is always suppressed. 
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0.5 



0.0 0.2 0.4 0.6 0.8 1.0 
n 

FIG. 14. The correlation exponent K p as a function of the bandfilling n for different values of 
U (U/t = 1, 2, 4, 8, 16 for the top to bottom curves). Note the rapid variation near n = 1 for small 
U. 



It should be emphasized here that the results of fig.[T4| are valid for n — > 1, but not 
for n = 1. In fact, in that latter case, there is a gap in the charge excitation spectrum, 
as expected from the umklapp term discussed below, and the correlations of (fi p become 
long ranged. Close to half-filling, the asymptotic behavior of the charge part of correlation 
functions like ( |4.42| ) is essentially determined by the motion of the holes. Writing the density 
of holes as p 



1/(1 — n) one then expects a crossover of the form ||1 06 
(n(x)n(O)) « cos(2A; F x)[l + {pxf] K " l2 x- 1 ln~ 3/2 (x) 



(6.40) 



for the 2k f part of the density correlation function, and similarly for other correlation 
functions. Clearly, only for x ^> l/p are the asymptotic power laws valid, whereas at 
intermediate distances 1 C i < 1//9 one has effectively K p = 0. Clearly, the form (|6.40 ) 
provides a smooth crossover as n — > 1. 

Results equivalent to the present ones can be obtained using the conformal invariance of 
the Hubbard model ||110| , |111|| . These results have subsequently be generalized to the case 
with an applied magnetic field ||112f |. The method for calculating the exponent K p is quite 



general and has been used for example for generalized Hubbard models ||105| , |113|| and the 
t-J model. rajTH 



3. Transport properties and the metal-insulator transition 

The exact solution of Lieb and Wu can also be combined with the long-wavelength 
effective Hamiltonian (|4.30| ) to obtain some information on the frequency-dependent con- 
ductivity a{uj). In particular, from eq. Q4.371) there is a delta function peak at zero frequency 



of weight <7 = 2K p u p , with the results plotted in fig.[L5[ As expected, far from half-filling, 
the dc weight is nearly independent of interaction strength, peak. On the other hand, for 
n — > 1 o"o vanishes linearly, implying a linear variation of the ratio n/m* with "doping", 
and for exactly half-filling the dc conductivity vanishes. There thus is a metal-insulator 
transition at n = 1, with the insulating state existing only at n = 1. One can also obtain 



the total frequency-integrated conductivity. [|105|| One then sees that as n — > 1 more and 
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0.0 0.2 0.4 0.6 0.8 1.0 
n 

FIG. 15. The weight of the dc peak in o~(iv) as a function of bandfilling for different values of 
U/t (U/t = 1,2,4,8,16 for the top to bottom curves). 

more weight is transferred from the zero-frequency peak to higher frequencies, essentially 
but not exclusively to frequencies above a "Mott-Hubbard gap" . 

It is interesting to discuss the metal-insulator transition in more detail. In the bosoniza- 
tion description, the insulating state at half filling is due to an extra term 

H u = J dxcos(V%) (6.41) 

which is due to umklapp scattering. The charge part of the Hamiltonian can then be 
transformed into a model of massive fermions, with energy-momentum relation = 
±(v 2 k 2 + A 2 ) 1 / 2 , ||116|| where A is the charge excitation gap created by umklapp scattering. 



At half-filling, this term leads to a gap A c in the charge excitation spectrum and therefore 
to insulating behavior. At half-filling all negative energy states are filled, all positive energy 
states are empty, and one has an insulator. Doping with a concentration n* of holes (n < 1), 
some of the negative energy states become empty and only states with \k\ > k* F oc n* are 
filled. Because of the vanishing interaction, a standard formula can be used ||117|| and gives 



a positive thermopower, i.e. approaching the metal-insulator transition from n < 1, the 
thermopower is hole-like, [ 118 , 119(1 whereas obviously far from the transition (n <C 1) it is 



electron-like. The exactly opposite behavior can be found for n > 1. 

The fact that u p and <To vanish linearly as n — > 1 seems to be consistent with a divergent 
effective mass at constant carrier density because u p ~ 1/m*, ctq « n/m*. A constant carrier 
density is also consistent with the fact that kp = im/2 is independent of U. It is not 
consistent with the hole-like sign of the thermopower as n — > 1 from below, nor with the 
electron-like sign as n — > 1 from above: if the carriers are holes, the carrier density is the 
density of holes: n* = 1 — n. Treating the holes as spinless fermions, as already mentioned 
before, one expects ctq — > because n* — > 0, and 7 — > 00 because the density of states 
of a one-dimensional band diverges at the band edges. This agrees with what was found 
explicitly in section |V1 B 2j. What is not so easily understood in this picture is the fact 



that kp (i.e. the location of the singularity of is given by its free-electron value 7m/2, 
rather then being proportional to n*. One should however notice that rtk is given by the 
single-particle Green's function, which contains both charge and spin degrees of freedom. 
The location of kp then may possibly be explained by phase shifts due to holon-spinon 
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interaction, 
solution 



This is in fact suggested by the structure of the wavefunction of the exact 
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The magnetic properties do not agree with what one expects from an effective mass 
diverging as n — > 1: u a and therefor \ remain finite. Moreover, the NMR relaxation rate 
would have the behavior l/T% = aT + (3yT, where the first (Korringa) term comes from 
fluctuations with q w 0, whereas the second term comes from antiferromagnetic fluctuations 
with q fa 2kp- None of these properties is strongly influenced by the diverging effective mass 
observed e.g. in the specific heat. This fact is of course a manifestation of the separation 
between spin and charge degrees of freedom. 

The discussion can be generalized to metal-insulator transitions near other rational band- 
fillings, i.e. one-third or one-quarter filled bands. [120| In the Hubbard model, however, an 
insulating phase occurs only at half-filling. 

It is worthwhile to point out here that the behavior of the one-dimensional system is 
quite different from that of scenarios for the metal-insulator transition in Fermi liquids. For 
example, in the "nearly-localized" description, [121] effective mass effects dominate, giving 



rise to a simultaneous divergence of specific heat and spin susceptibility. On the other hand, 
in infinite dimension |93|,|94| a behavior similar to the one-dimensional case is observed: 
approaching the insulating state by varying the bandfilling the spin susceptibility remains 
finite but the specific heat coefficient diverges. ||122| 



4- Spin-charge separation 

The Hubbard model also provides a rather straightforward interpretation of the spin- 
charge separation discussed above. Consider a piece of a Hubbard chain with a half-filled 
band. Then for strong U there will be no doubly-occupied sites, and because of the strong 
short-range antiferromagnetic order the typical local configuration will be 

■■■nnnnnn--- 

Introducing a hole will lead to 

•■■ tint o Tint! •■■ 

and after moving the hole one has (note that the kinetic term in the Hamiltonian does not 
flip spins) 

-uoumuu- 

Now the hole is surrounded by one up and one down spin, whereas somewhere else there are 
two adjacent up spins. Finally, a few exchange spin processes lead to 

-uouuum- 

Note that the original configuration, a hole surrounded by two up spins has split into a 
hole surrounded by antiferromagnetically aligned spins ("holon") and a domain-wall like 
configuration, two adjacent up spins, which contain an excess spin 1/2 with respect to the 
initial antiferromagnet ("spinon"). The exact solution by Lieb and Wu contains two types 
of quantum numbers which can be associated with the dynamics of the spinons and holons, 
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respectively. We thus notice that spinons and holons ||123| , |124|| have a well-defined meaning 
in the present one-dimensional case. 

The above pictures suggest that, as far as charge motion is concerned, the Hubbard 
model away from half-filling can be considered as a one-dimensional harmonic solid, the 
motion of the holes providing for an effective elastic coupling between adjacent electrons. 
This picture has been shown to lead to the correct long-distance correlation functions for 
spinless fermions P7|jl25| |. For the case with spin, this suggests that one can consider the 
system as a harmonic solid with a spin at each site of the elastic lattice (lattice site = 
electron in this picture), and this gives indeed the correct spin correlations 



VII. CONCLUSION AND OUTLOOK 



In these notes, the different behavior of interacting many-fermion systems in one and 
higher dimension have been contrasted. In dimension bigger than one, Landau's Fermi liquid 
theory is based on the postulated existence of quasiparticles which have properties closely 
related to those of a noninteracting system. On the other hand in one dimension, essentially 
because of the geometric constraints on the particle dynamics no well-defined quasiparticle 
excitations exist. Instead, the low-lying excitations are collective spin or charge density 
oscillations, giving rise to the interesting phenomenon of spin-charge separation, as well as 
to non-universal power law behavior of correlation functions. This Luttinger liquid physics 
is most appropriately described by the bosonization formalism developed in sec.[IV[ We may 
note here that bosonization can equally be applied in higher dimension, ||126| , |127|| however 
at least for "normal" interactions one finds Fermi liquid properties. 

It is clearly of interest to investigate the possibility of a crossover between Luttinger and 
Fermi liquid. A most natural candidate for such behavior are systems of parallel chains 
coupled by a small interchain hooping integral t±. Such a model is appropriate for the 
discussion or real gimsz-one-dimensional conductors. This type of models has been studied 
for some time but there is still a number of open questions. Nevertheless, most work points to 
a crossover between Luttinger (at high temperature) and Fermi liquid (at low temperature) 
behavior at a temperature determined by t±. [|128| , |1 29| . |]~^ , pL 3Q|JT~3Ty] At lower temperatures 
three-dimensionally ordered phases can appear. A different scenario has been advocated 
by Anderson and collaborators ||132| , |133[ | who argue that even well below T ~ t± particle 
movement in the direction transverse to the chains remains incoherent, i.e. there would still 
be no three-dimensional Fermi liquid. Another approach is to consider the spatial dimension 
as a continuously varying parameter. ||134|| One then finds that the Luttinger liquid disappear 
as soon as d > 1. 

It is worthwhile to mention here that the question of Fermi liquid or non-Fermi liquid 
is of interest in a number of other contexts. For example, the low-temperature properties 
of a "normal" Kondo impurity can be described as a local Fermi liquid. ||135|| On the other 
hand, the so-called "overscreened" Kondo effect ||136|| is characterized by non-Fermi liquid 
behavior, in particular a divergent local spin susceptibility. ||137|| A review of this subject 
has been given by Emery and Kivelson. ||13£ 



The other major area of current interest in possible non-Fermi liquid behavior is the 
physics of strongly correlated systems in two dimensions, and in particular of the high- 
T c superconducting copper oxides. Here a number of experimental results in the normal 
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conducting state suggest non-Fermi liquid behavior. [|1 39| |1 4 1|| The theoretical situation 
is much less clear. On the one hand, a number of radically non-Fermi liquid pictures 
have been proposed, e.g. a two-dimensional Luttinger liquid ||140|| or the so-called gauge 
theory approach. [ |142|| On the other hand, there are also models with only relatively minor 
modifications of the Fermi liquid picture. A rather complete overview of the rather unclear 
current situation (1994) can be found in recent conference proceedings. ||143|| Finally, an 
interesting non-Fermi liquid state has been proposed for the half-filled Landau level in two 
dimensions. 1144 
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